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EXECUTIVE  PROGRAM  SUMMARY 


This  NRL  Memorandum  Report  describes  the  research  accomplishments  of 
the  Radiation  Hydrodynamics  Branch,  in  part  in  collaboration  with  Sandia 
National  Laboratories,  during  Fiscal  Year  1991.  It  consists  of  6  sections 
whose  authors  are  separately  identified  at  the  beginning  of  each  section. 

Perhaps  the  most  noteworthy  achievement  this  past  year  was  the  first- 
ever  demonstration  of  a  population  inversion  in  the  x-ray  region  pumped 
solely  by  resonant  photoexcitation.  Last  year's  report  discussed  a  time  - 
integrated  neon  spectrum  which  suggested  inversion.  This  year's  first 
section  presents  detailed  time-  resolved  spectra  of  He  -  like  neon  pumped  by 
both  (resonant)  sodium  and  (non-resonant)  magnesium,  which  spectra  prove  the 
existence  of  the  inversion.  Also,  this  achievement  (using  Sandia  National 
Laboratories'  10  -  megampere  Saturn  accelerator)  represents  the  first 
population  inversion  driven  by  a  Z-pinch.  The  major  goal  of  the  continuing 

O 

Saturn  experiments  is  to  detect  and  measure  the  gain  of  the  231  A  laser  line. 
We  believe  the  feasibility  of  creating  x-ray  inversions  by  resonant 
photoexcitation  has  now  been  demonstrated.  Success  in  demonstrating  gain 
will  open  the  door  to  a  new  and  radically  different  class  of  x-ray  lasers 
with  large  output  energies.  As  this  report  was  being  prepared,  we  received 
the  news  that  the  material  in  Sec.  I  has  been  accepted  for  publication  in  the 
Feb.  10,  1992  Physical  Review  Letters. 

The  following  two  sections  of  this  report  are  also  concerned  with 
resonantly  photopumped  ultrashort  wavelength  lasers,  but  of  considerably 
different  design.  Sec.  II  presents  the  first  detailed  analysis  of  the  A1  XI 
-  Mg  IX  scheme  proposed  by  Krishan  and  Trebes  in  1984.  Although  the  atomic 
physics  of  these  ions  is  more  complicated  than  that  of  the  simple  He  -  like 
stages  of  sodium  and  neon,  we  conclude  that  this  scheme  is  as  promising  as 
Na-Ne,  and  might  well  be  suitable  as  a  "tabletop"  device  for  smaller  machines 
because  of  the  ease  of  fabrication  of  Al  -  Mg  targets.  In  Sec.  Ill,  a  class 
of  lasers  featuring  "automatic"  pump  resonance  is  discussed.  In  this  method, 
a  hot,  dense  neonlike  krypton  shell  pumps  a  cooler,  more  tenuous  lasant  inner 
krypton  plasma.  The  advantags  of  the  scheme  are  that  all  the  pump  lines 
automatically  match  in  wavelength  (no  accidental  coincidence  is  required  as 
for  Na-Ne  or  Al  -Mg),  and  ,  the  low  density  interior  plasma  is  configured  to 
minimize  x-ray  refraction  effects  which  degrade  beam  quality.  Large  output 
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energies  and  long  lasing  times  (several  nsec  or  more)  are  additional 
attractive  features  of  this  laser  when  configured  on  an  accelerator  as  a 
coaxial  Z  pinch. 

One  of  the  most  exciting  new  technologies  which  has  emerged  in  the  past 
3  years  or  so  is  that  of  the  high-pover,  femtosecond-duration  optical  laser. 
Potential  applications  are  numerous  and  have  been  widely  discussed  in  the 
literature.  The  current  conventional  wisdom  is  that  such  devices  are  the 
number  one  candidate  for  producing  cheap,  room-sized  "tabletop"  x-ray 
lasers.  Section  IV  presents  a  detailed  study  of  radiation  energetics  and 
target  coupling  for  the  interaction  of  an  ultrabright  laser  pulse  with  a 
planar  fluorine  target.  Section  V  describes  the  development  and  validation 
of  the  first  fully  time-dependent  radiative  transfer  model  which  can  be 
employed  in  analysis  of  laser  plasmas.  Since  the  velocity  of  light  is  0.3  ym 
/  fsec,  in  many  of  the  ultrashort  pulse  laser  heated  plasmas,  the  light 
travel  time  is  comparable  to  the  heating  or  cooling  time.  It  is  anticipated 
that  this  model  will  be  coupled  to  a  realistic  ,  multilevel  atomic  data  base 
in  FY92.  Taken  together,  the  capabilities  detailed  in  Secs  IV  and  V 
constitute  a  firm  analytical  base  for  the  design  of  an  eventual  compact, 
high-output  x-ray  laser. 

In  the  concluding  Sec.  VI,  one  of  the  most  baffling  continuing 
mysteries  of  x-ray  laser  experiments  is  rigorously  analyzed.  Ever  since 
neonlike  x-ray  lasers  were  first  demonstrated  in  1984,  it  has  been  generally 
observed  that  the  strongest  predicted  laser  transition,  that  arising  from  an 
excited  J=0  level,  is  anomalously  weak.  Many  explanations  have  been  offered, 
with  as  yet  no  consensus  developed  in  the  community.  Sec  VI  presents  a  near- 
exhaustive  analysis  of  the  possibilities  which  are  not  yet  scientifically 
ruled  out.  While  this  major  contribution  cannot  be  adequately  summarized 
here,  one  of  the  authors'  most  important  findings  concerns  the  effect  of 
inability  to  calculate  collision  rate  coefficients  to  better  than  factors  of 
two  or  three.  The  authors  have  discovered  several  instances  in  the 
literature  where  the  accuracy  of  a  key  monopole  collision  rate  has  been 
seriously  challenged,  contrary  to  many  existing  claims  of  excellent  accuracy 
for  this  fundamental  atomic  parameter.  Their  study,  among  other  approaches, 
varies  this  rate  through  plausible  values  and  concludes  that  this  rate  must 
be  measured,  and,  if  lover  than  the  present  consensus  values,  could  explain 
this  persistent  mystery. 
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Abstract 

The  broadband  radiation  emitted  from  a  Na  z  pinch  is  used  to  photoionize  Ne 
to  the  He-like  ground  state  and  radiation  from  the  Na  1s2-1s2p  ^  Pi  transition  is 
used  to  resonantly  photoexdte  the  Ne  1s2-1s4p  ^P-\  transition.  Time-resolved 
and  time-integrated  spectral  measurements  of  the  Ne  K-shell  emission 
demonstrate  the  first  population  inversion  driven  by  a  z  pinch.  This  is  the  first 
experiment  in  any  medium  to  demonstrate  a  soft  x-ray  inversion  pumped  solely 
by  resonant  photoexcitation. 

Article 

Research  conducted  during  the  past  10  years  has  demonstrated  that  pulsed- 
power  drivers  are  capable  of  efficiently  creating  very  energetic  radiation 
sources  with  a  variety  of  spectral  distributions.^  Many  have  suggested  that 
these  drivers  can  be  used  for  applications  such  as  creating  photopumped  x-ray 
lasers^,  imploding  inertial  confinement  fusion  capsules,^  and  studying  the 
photoionization  kinetics  of  plasmas  in  intense  radiation  fields.^  A  serious 
potential  problem  with  using  pulsed-power  drivers  is  that  the  powerful  electrical 
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pulse  (currents  of  many  MA  and  voltages  of  many  MV)  flowing  within  millimeters 
of  the  various  "targets"  may  obscure  the  effects  of  the  radiation  drive  and 
hopelessly  complicate  the  execution  of  this  class  of  experiment.  In  this  Letter 
we  show  that  these  problems  are  separable  and  present  the  first  unambiguous 
experimental  demonstration  of  population  inversion  driven  by  a  z  pinch,  and  the 
first  soft  x-ray  inversion  pumped  solely  by  resonant  photoexcitation  in  any 
medium. 

We  are  performing  photopumping  experiments  using  the  sodium/neon 
resonant  photoexcitation  x-ray  laser  scheme.  This  scheme,  which  has  received 
extensive  theoretical  study,®"®  employs  radiation  from  the  1s2-ls2p  ^P-| 

transition  at  11.0027  A  in  He-like  Na  to  resonantly  photoexcite  the 
Ne  1s2.ls4p  ^Pi  transition  at  11.0003  A  in  Hs-like  Ne.  Electron  and  ion 

collisions  are  predicted  to  transfer  a  large  fraction  of  the  excited  4p  singlet 
population  to  the  4d  and  4f  levels.  The  line  calculated  to  have  the  highest  gain 
is  the  1s3d  ^D2  -  1s4f  transition  at  231  A.  This  resonance  is  attractive 

because  of  the  excellent  match  between  the  Na  and  Ne  resonant  transitions, 
the  intrinsic  strength  of  the  Na  pump  line,  and  the  relative  ease  in  creating  a 
large  population  of  Ne  ions  in  the  He-like  ground  state. 

To  produce  a  bright  Na  pump  line  one  needs  a  high  density 
(>  10^®  ions  cm"®),  high  temperature  (>  300  eV)  Na  plasma  in  order  to  both 
maximize  the  radiation  rates  and  insure  that  a  large  population  of  Na  ions  are  in 
excited  He-like  states.  The  Ne  lasant  plasma,  however,  needs  to  be  at  a  lower 
density  (*10^®  ions  cm"®)  and  temperature  (<  100  eV)  in  order  to  both 
maximize  the  population  of  Ne  ions  in  the  He-like  ground  state  and  minimize 
collisional  excitation  of  the  lower  lasing  level. 

In  our  experimental  arrangement  (see  Fig.  1)  such  divergent  conditions  have 
been  achieved.  We  create  the  Na  plasma  by  imploding  an  array  of  16,  20-mm- 
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long,  75-p.m-diameter  Na  wires  using  the  Saturn  accelerator^  as  a  z-pinch 
driver.  Saturn  implodes  the  Na  wire-array  load  with  a  peak  current  of  10  MA  in 
an  electrical  pulse  with  a  peak  power  of  20  TW  and  a  40-nsec  FWHM.  The  total 
broadband  radiation  output  from  this  plasma  is  measured  to  be  400  kJ  in  a 
pulse  with  a  peak  power  of  10  TW  and  a  40-nsec  FWHM.  Measurements  show 
that  90%  of  this  radiation  is  emitted  at  energies  less  than  1  keV.  Since  the  final 
pinch  diameter  is  observed  to  be  2.5  mm  (implying  ion  densities  of 
4  X  10^®  cm'3),  we  assume  that  the  broadband  radiation  from  this  plasma 
can  be  approximated  by  the  emission  from  a  89-eV  blackbody  (simply  the 
blackbody  equivalent  temperature  for  the  measured  radiating  area  and  power). 
6  kJ  of  energy  is  typically  observed  to  be  radiated  in  the  1 1  A  Na  pump  line  in  a 
pulse  with  a  peak  power  of  200  GW  and  a  20-nsec  FWHM.  This  pulse  is 
basically  Gaussian  in  shape  until  about  25  nsec  after  the  peak,  at  which  time  a 
"foot"  appears  that  slowly  decays  away.  Approximately  16%  of  the  energy  and 
power  which  is  radiated  between  1-3  keV  appears  in  the  11  A  pump  line.  The 
pump  line-width  has  been  measured^  ^  with  a  high-resolution  time-integrating 
spectrometer  and  found  to  be  approximately  20  mA.  Calculations  of  this  opacity 
broadened  line  width,  based  on  the  above-cited  diameter  and  density  and  the 
calculated  range  of  He-like  Na  fraction,  result  in  a  predicted  line  width  of 
10-20  mA.  Since  the  ionization  balance  and  thus  the  opacity  of  this  line  vary 
with  time,  a  definitive  measurement  of  the  pump  line  width  will  require  a  high- 
resolution  time-resolving  spectrometer.  This  difficult  measurement  has  not 
been  made  to  date. 

We  use  the  broadband  radiation  from  the  Na  pinch  to  photoionize  the  Ne  to 
the  He-like  ionization  state.  The  ionization  balance  in  the  Ne  is  dominated  by 
the  Na  z-pinch  radiation  field  and  not  the  Ne  electron  temperature. 
Independent  simulations  of  these  experiments  by  Sandia  National 
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Laboratories,  the  Naval  Research  Laboratory,^  ^  and  Lawrence  Livermore 
National  Laboratory'*  ^  show  that  53-80%  of  the  Ne  ions  are  in  the  He-like 
ground  state  even  though  the  Ne  electron  temperature  is  only  25-30  eV. 

The  target  used  in  these  experiments  is  located  20  mm  away  from  the 
z-pinch  axis  (outside  the  current  return  posts)  and  consists  of  Ne  at  a  pressure 
of  10  Torr  (4  x  10^^  atoms  cm‘3)  contained  within  a  "gas  cell"  target.  The 

viewing  and  illumination  windows  of  this  target  are  made  of  free  standing, 
5000-A-thick  Lexan  (C10H14O3).  The  illumination  window  which  faces  the  Na 

pinch  is  10-mm  long  and  3-mm  wide.  The  two  viewing  windows  which  face  the 
lasant  spectrometers  are  circular  in  shape  with  a  diameter  of  6  mm.  Radiation- 
hydrodynamic  calculations  predict  that  the  Lexan  illumination  window  is  heated 
by  the  broadband  radiation  from  the  Na  pinch  and  expands,  compressing  a 
few-hundred-micron-thick  region  of  the  Ne  to  densities  suitable  for  lasing 
(*10^®  ions  cm”3)  at  the  time  of  peak  Na  11  A  pump  emission.  This 
compressed  region  is  approximately  one  mean  free  path  thick  to  the  1 1  A  pump 
line.  Not  only  is  this  target  both  designed  and  measured  to  achieve  suitable 
laser  conditions,  but  as  discussed  later,  our  measurements  indicate  that  it 
contains  no  current  and  associated  instabilities.  This  eliminates  a  major 
difficulty  of  pulsed-power  driven  x-ray  lasers.  While  the  present  20  mm 
separation  facilitates  clean  diagnosis  and  is  sufficient  to  demonstrate 
population  inversion,  the  small  solid  angle  subtended  by  the  sodium  flashlamp 
can  be  significantly  improved  in  future  experiments.  The  use  of  more  advanced 
targets  could  lead  to  gain  in  the  4-2  and  3-2  transitions  which  have  shorter 
wavelengths  and  higher  quantum  efficiency. 

The  total  broadband  radiated  energy  and  power  from  the  Na  pinch  is 
measured  with  an  unfiltered,  thin-film  resistive  bolometer.  The  radiated  energy 
and  power  above  1  keV  is  measured  with  a  filtered  bolometer,  a  time- 
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integrating  convex  KAP  crystal  spectrometer,  and  an  array  of  three  absolutely- 
calibrated.  filtered  x-ray  diodes A  12-frame,  time-resolving,  x-ray  pinhole 
camera  is  used  to  monitor  the  implosion  quality  of  the  Na  pinches. 

The  K-shell  emission  between  9-14  A  from  the  Ne  lasant  is  monitored 
simultaneously  with  time-resolving  and  time-integrating  elliptical  crystal  (RAP) 
spectrometers  (see  Fig.  1).  The  time-resolving  spectrometer  is  located  1.2  m 
from  the  target  and  uses  a  7-frame,  microchannel-plate-intensified.  x-ray 
framing  camera  as  the  detector.  The  time-integrating  spectrometer  is  located 
0.35  m  from  the  target  and  uses  Kodak  2497  x-ray  film.  We  insert  metal  tubes 
between  the  x-ray  laser  target  and  each  spectrometer  to  guard  against  any  Na 
contaminant  radiation  entering  the  iine-of-sight  of  these  spectrometers.  The 
ends  of  these  baffles  are  placed  inside  counterbored  holes  located  at  both  the 
x-ray  laser  target  and  the  entrance  port  of  the  spectrometers.  To  insure  that  no 
Na  contaminant  radiation  was  observed  by  these  lasant  spectrometers,  we 
conducted  a  "null"  shot  by  replacing  the  Ne  with  N  at  10  Torr.  A  nominal  Na 
pinch  was  produced  on  this  shot  and  no  fines  were  observed  on  either  of  the 
lasant  spectrometers. 

We  looked  for  an  enhancement  of  the  He-like  Ne  1s^-1s4p  line  relative  to 
the  He-like  Ne  1s^-1s3p  line  to  demonstrate  that  the  Na  pinch  photopumps  the 
Ne  plasma.  The  three  independent  kinetics  calculations  have  shown  that  if  an 
inversion  exists  between  these  two  transitions  (i.e.  the  1s-4p  line  is  brighter  than 
the  1s-3p  line),  then  at  the  Ne  densities  encountered  in  these  experiments,  an 
inversion  also  exists  between  the  1s4f  and  1s3d  singlet  levels.  The  exact 
magnitude  of  the  gain,  however,  is  uncertain.  The  kinetics  calculations  differ  in 
the  predicted  He-like  fraction,  oscillator  strengths  of  the  pumped  and  laser 
transitions,  width  of  the  laser  transition,  and  level  of  detail  in  modeling  the 
excited  states.  Uncertainties  of  10-60%  in  these  factors  lead  to  the  very 
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sensitive  prediction  of  gain  varying  from  0.2  to  1.5  cm''*.  Additional 
experimental  uncertainties  arise  from  the  width  of  the  11  A  pump  line  at  peak 
power  and  the  details  of  the  soft  photoionizing  spectrum. 

To  ensure  that  the  enhanced  emission  observed  in  the  1s^-1s4p  line  is  due 
to  resonant  photopumping  by  the  Na  11  A  pump  line,  we  performed  a  second 
type  of  "null"  shot.  In  this  "null"  shot  the  Na  wire-array  load  was  replaced  with  a 
Mg  wire-array  load  in  an  attempt  to  generate  the  same  broadband  and  general 
keV  radiation  as  a  Na  pinch  but  eliminate  the  1 1  A  Na  pump  line. 

Figure  2a  shows  typical  time-integrated  and  time-resolved  Ne  spectra  from  a 
shot  with  a  Na  pinch  and  a  Ne-filled  target.  The  time-integrated  spectrum  has 
been  corrected  for  the  spectrometer  response^  ^  (using  tabulated  values  for  the 
film  response,  filter  transmission,  crystal  reflectivity,  and  spectrometer  geometry) 
and  background  light.  In  addition,  the  time-resolved  spectra  were  also 
corrected  for  the  gold  photccathode  response  and  the  variation  in  gating 
voltage  (and  hence  detector  gain)  between  frames.  For  this  shot,  each  gating 
pulse  was  4  nsec  wide  and  spaced  5  nsec  apart  (center  to  center).  Since  the 
framing  camera  was  run  at  very  high  gain  (=1  kV  across  the  microchannel 
plate),  the  bright  1s2-ls2p  line  might  be  slightly  saturated  on  some  frames.  The 
high  gain  also  resulted  in  "noisier"  time-resolved  data  than  time-integrated  data. 

As  can  be  seen  from  Fig.  2a,  the  He-like  Ne  1s2-1s4p  line  is  about  twice  as 
bright  as  the  Is^-lsSp  line,  demonstrating  an  inversion  between  the  1s4p  and 
1s3p  singlet  levels.  Integration  of  these  spectra  gives  an  intensity  ratio  for  the 
1s2-1s4p  line  to  the  1s^-1s3p  line  of  1.7  for  the  time-integrated  spectrum  and  a 
ratio  of  2.4  at  the  time  of  peak  emission  of  the  Na  1 1  A  pump  line.  This  pattern 
was  observed  on  each  of  three  Na-pumped  shots.  The  three  independent 
simulations  predict  this  line  ratio  to  be  between  1.4  and  2.2  at  the  time  of  peak 
Na  11  A  pump  emission.  The  spontaneous  decay  rate  of  the  1s^-1s3p  line 
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exceeds  that  of  the  1s2-is4p  transition  by  a  factor  of  2.4.  Even  with  the  effects 
of  self-absorption,  population  inversion  is  required  to  make  the  4p  transition 
twice  as  bright  as  the  3p,  as  is  observed  here. 

The  excited  state  density  of  the  pumped  1s4p  level,  N(4p).  may  be  inferred 
from  the  measured  power  P  (photons  sr'^  sec'^)  of  the  11  A  Ne  line  from 
N(4p)=4nP/AV  where  A  is  the  line's  decay  rate  (10^^  sec"^)  and  V  is  the 
observed  Ne  volume.  By  normalizing  the  absolute  energy  measured  with  the 
time-integrating  spectrometer  to  the  relative  time  history  measured  with  the 
time-resolving  spectrometer,  we  estimate  that  P=4  x  10^®  photons  sr'^  sec'^. 
The  observed  volume  is  the  product  of  the  observed  window  length  (3  mm),  the 
depth  of  penetration  of  the  Na  pump  radiation  into  the  Ne,  and  the  depth  into  the 
Ne  seen  by  the  time-integrating  instrument.  The  latter  two  quantities  are 
inversely  proportional  to  the  line  optical  depth  and  therefore  to  the  He-like 
fraction  f.  As  a  result,  N(4p)  is  proportional  to  f^.  if  f=0.53  is  assumed  (LLNL 
model.  Ref.  12),  we  infer  N(4p)=2.0  x  10^^  cm*^,  compared  to  the  Ref.  12 
calculation  of  3.0  x  10^^  cm*^.  However,  the  f  of  0.8  predicted  by  the  NRL  and 
SNL  models  leads  to  N(4p)=4.5  x  10^^  cm'^,  compared  to  the  corresponding 
model  calculations  of  5.9  x  10^^  and  2.9  x  10^^  cm’^.  Therefore,  taking  into 
account  the  dependence  of  the  viewing  volume  on  t,  the  measured  power  of  the 
photopumped  line  is  consistent  with  all  3  models  to  within  35%. 

Another  interesting  feature  of  the  time-resolved  spectra  shown  in  Fig.  2a  is 
the  general  time  history  of  the  Ne  K-shell  lines.  Note  that  no  lines  appear  until 
frame  3  of  the  time-resolved  data.  If  the  ionization  balance  and  the  excitation  of 
the  Ne  plasma  were  dominated  by  stray  currents  or  by  high  energy  electron 
beams  flowing  in  this  plasma,  then  one  would  expect  these  lines  to  appear 
much  earlier  in  time  since  this  shot  had  a  65  nsec  implosion  time  (time  from 
10%  of  maximum  current  to  stagnation).  The  fact  that  these  K-shell  lines  follow 
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the  pinch's  keV  radiation  pulse  (and  not  the  current  pulse)  indicates  that  both 
the  ionization  balance  and  the  excitation  of  the  Ne  plasma  are  dominated  by  the 
radiation  field  produced  by  the  Na  pinch. 

Figure  2b  shows  the  time-integrated  and  time-resolved  Ne  spectra  from  the 
Mg/Ne  "null"  shot.  For  this  shot  the  gating  pulses  on  the  time-resolving 
spectrometer  were  8  nsec  wide  and  spaced  5  nsec  apart  (center  to  center), 
resulting  in  some  overlap  between  adjacent  frames.  The  main  feature  to  note 
from  this  spectrum  is  the  fact  that  the  Is^-lsAp  line  is  no  longer  brighter  than  the 
Is^-isSp  line.  The  spectrally  integrated  intensity  ratio  of  the  1s2-1s4p  line  to 
the  ts^-lsSp  line  is  1.0  for  both  the  time-integrated  spectrum  and  frame  number 
4  of  the  time-resolved  spectra.  This  observation,  combined  with  the  fact  that  no 
Na  radiation  was  observed  on  the  Na/N  "null"  shot,  confirms  our  interpretation  of 
resonant  photopumping  as  the  cause  of  the  observed  inversion  between  the 
1s2.is4p  and  Is^-lsSp  lines  in  the  Na/Ne  shots.  Another  interesting  feature  of 
the  Ne  spectra  from  the  Mg/Ne  shot  is  the  relative  brightness  of  the  Is^-lsSp, 
1s2-1s4p,  and  Is^-lsSp  lines.  The  integrated  line  intensities  of  these  spectral 
lines  are  in  the  ratio  1.0:1 .7:1 .7  for  the  time-integrated  spectrum  and  1.0:1 .3:1.3 
for  frame  number  4  of  the  time-resolved  spectra.  This  is  suggestive  of 
recombination  from  H-like  Ne  ions  in  the  cold,  moderate  density  Ne  plasma.  A 
population  of  H-like  ions  is  known  to  exist  on  this  shot  since  the  H-like  Ne  1s-2p 
resonant  line  is  observed  on  frames  3-6  of  the  time-resolved  spectra. 

It  is  Important  to  note  that  no  H-like  lines  have  ever  been  observed  on 
Na-pumped  shots.  Moreover,  the  1s^-1s5p  line  is  absent  when  a  Na  plasma  is 
the  pump.  It  can  be  shown  that  the  three-body  recombination  rate  scales  as 
n5.33  for  the  plasma  conditions  discussed  here,  where  n  is  the  principal 
quantum  number.  This  scaling  is  deduced  from  a  detailed  balancing  of 
collisional  ionization  (using  the  Lotz  ionization  formula^®)  with  three-body 
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recombination.  The  absence  of  n=5  radiation  during  the  Na*pumped  shots,  as 
well  as  the  lack  of  H-like  lines  or  any  recombination  continuum,  rules  out 
recombination  as  a  pump  mechanism  for  the  inversion.  In  contrast,  the 
Mg-pumped  shots  exhibit  the  H-like  Ne  1s-2p  line  and  the  recombination 
continuum,  as  well  as  a  strong  He-like  Ne  1s^-1s5p  line,  which  are  all 
characteristic  of  recombination.  The  keV  lines  from  Mg  are  all  of  sufficient 
energy  to  ionize  He-like  Ne  and  were  also  emitted  with  40%  more  power  than 
those  of  Na,  fully  consistent  with  the  above  picture.  Analysis  of  the  (somewhat 
noisy)  slope  of  the  recombination  continuum  of  the  Mg  shot  yields  an  electron 
temperature  of  30  eV,  consistent  with  the  radiation-hydrodynamic  simulations  of 
these  experiments."’"’ 

In  summary,  we  have  made  time-integrated  and  time-resolved  spectral 
measurements  which  demonstrate  population  inversion  by  resonant 
photopumping  in  the  Na/Ne  x-ray  laser  system.  Success  in  demonstrating  gain 
on  this  or  other  laser  transitions  would  open  the  door  to  a  class  of  x-ray  lasers 
with  large  output  energies  and  lasing  times  of  several  nsec  or  longer. 

We  would  like  to  thank  B.  Hammel  for  the  many  early  discussions  which 
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Saturn  operations  crew.  The  targets  were  supplied  by  Luxel  Inc.  of  Friday 
Harbor,  WA.  The  Sandia  authors  would  also  like  to  thank  J.  P.  VanDevender  for 
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SDIO/T/IS  and  the  U.S.  Department  of  Energy  under  Contract 
DE-AC04-76DP00789. 


9 


References 


^  N.  R.  Pereira  and  J.  Davis,  J.  Appl.  Phys.  64,  R1  (1988). 

^  M.  K.  Matzen  etal.,  J.  Phys.  Colloq.  (Paris)  47,  C6-135  (1986). 

^  J.  P.  VanDevender  and  D.  L  Cook,  Science  232,  831  (1986). 

^  M.  J.  Eckart,  private  communication. 

®  A.  V.  Vinogradov,  1. 1.  Sobelman,  and  E.  A.  Yukov,  Sov.  J.  Quantum 
Electron.  5,  59  (1975). 

®  P.  L  Hagelstein,  Plasma  Phys.  25, 1345  (1983). 

^  J.  P.  Apruzese  and  J.  Davis,  Phys.  Rev.  A.  31,  2976  (1985). 

®  S.  J.  Stephanakis  etal.,  IEEE  Trans.  Plasma  Sci.  16,  472  (1988). 

®  R.  B.  Spielman  etal..  1989  Proc.  2nd  Int.  Conf.  on  Dense  Z  Pinches,  eds.  N. 
R.  Pereira,  J.  Davis,  and  N.  Rostoker  (New  York,  American  Institute  of  Physics) 
pp  3-16. 

P.  Burkhalter  etal.,  in  Record  1991  IEEE  Int.  Conf.  Plasma  Sci., 
Williamsburg,  VA  ,1991,  p.  212. 

^  ^  J.  P.  Apruzese  et  a!.,  in  Proceedings  of  the  2nd  International  Colloquium 
on  X-Ray  Lasers,  York,  England,  ed.  G.  J.  Tallents  (Bristol,  lOP)  pp.  39-42 
(1991). 

J.  Nilsen  and  E.  Chandler.  Phys.  Rev.  A  44,  4591  (1991). 

F.  C.  Young,  S.  J.  Stephanakis,  and  V.  E.  Scherrer,  Rev.  Sci.  Instrum.  57, 
2174  (1986). 

B.  A.  Hammel  and  L.  E.  Ruggles,  Rev.  Sci.  Instrum.  59,  1828  (1988). 

B.  L.  Henke  and  P.  A.  Jaanimagi,  Rev.  Sci.  Instrum.  56.  1537  (1985). 

■>6  w.  Lotz,  Z.  Phys.,  216,  241  (1968). 


10 


ACHIEVABLE  PUMP  POWER  AND  GAIN  IN  THE  Al  XI  -  Mg  IX 
PHOTORESONANT  X-RAY  LASER 


J.  P.  Apruzese  and  M. 


Buie* 


Radiation  Hydrodynamics  Branch 
Plasma  Physics  Division 
Naval  Research  Laboratory 
Washington,  D.  C.  20375 


Abstract  -  The  aluminum-magnesium  photoresonant  x-ray  laser  scheme  is 

investigated  using  detailed  atomic  models  for  both  elements.  The 

aluminum  pump  source  is  assumed  to  be  a  Z  pinch  driven  by  a  current 

of  several  MA.  It  is  found  that  monochromatic  pump  powers  of  100 

0 

GW/eV  are  obtainable  within  the  48  A  aluminum  pump  line  over  a  broad 
range  (125-400  eV)  of  pump  plasma  temperatures.  The  gain  in  Mg  IX  is 
calculated  assuming  separations  of  1  or  2  cm  from  the  pump  plasma. 
The  highest  gain  appears  in  the  2s3d  ■^D2  -  2s4f  line  at  228  A. 

However,  photoionization  as  well  as  resonant  photoexcitation  appears 
to  be  crucial  to  the  success  of  this  scheme.  If  photoionization  by 
lines  other  than  the  pump  line  is  neglected,  Mg  IX  does  not  exist 
below  -  60  eV,  In  the  60-100  eV  Mg  plasma  temperature  range, 
collisional  excitation  from  low-lying  states  to  the  lower  laser  level 
severely  degrades  the  gain.  However,  when  photoionization  of  Mg  by 
all  of  the  strong  Al  XI  lines  is  taken  into  account.  Mg  IX  is  found 
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to  exist  at  much  lower  temperatures,  leading  to  respectable  gain  of 
3  cm”^  at  1  cm  Al-Mg  separation. 

*  Current  address:  Dept,  of  Nuclear  Engineering,  Univ.  of  Michigan, 

Ann  Arbor,  Michigan  48109 
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I. 


INTRODUCTION 


In  1984  Krishnan  and  Trebes  analyzed^  a  series  of  six  Be-like  ions  whose 
2s'‘  -  2s4p  transitions  are  closely  resonant  with  various  possible  pump 
lines  in  ions  ranging  from  Mn  VI  to  Al  XI.  The  Be-like  lasant  ions  range  in 
atomic  number  from  6  (C  III)  to  12  (Mg  IX).  By  exposing  a  plasma  containing  the 
Be-like  ion  to  the  resonant  radiation  from  the  pump  ion,  gain  in  4-3  singlet 
transitions  may  be  achieved  through  pumping  of  the  ls^2s4p^Pj^  level.  The 
lowest-atomic-number  resonant  pair,  C  III  -  Mn  VI  has  already  been  successfully 
demonstrated'^,  resulting  in  the  shortest  laser  wavelength  (2163  A)  yet  achieved 
by  resonant  photoexcitation,  with  the  possible  exception  of  the  recent  neonlike 
titanium  experiments  at  the  University  of  Rochester  and  Lawrence  Livermore 
National  Laboratory^.  The  highest-atomic-number  resonant  pair,  Al  XI-Mg  IX,  had 
been  previously  proposed  by  Seely^  and  discussed  by  Dixon  and  Elton. ^ 

Of  course,  many  other  resonant  photoexcitation  schemes  to  obtain  gain  in 
the  soft  x-ray  region  have  been  proposed^"^^.  One  of  the  most  attractive  and 
highly  developed  of  these  is  the  11  A  resonance  '  ^  between  the  Na  X  Is  -ls2p 
^Pj  and  Ne  IX  Is^  -  ls4p  ^Pj  lines  where  experimental  evidence  of  both 
fluorescence  and  population  inversion  has  been  produced^^’ Previous 
analyses  of  this  sodium-neon  scheme  have  emphasized  the  need  for  considerably 
different  conditions  in  the  pumping  sodium  and  pumped  neon  plasmas.  Indeed, 
conditions  in  the  experiments  demonstrating  fluorescence  and  population 
inversion  are  consistent  with  the  requirement  for  a  sodium  plasma  much  hotter 
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and  denser  than  the  neon.  The  Be-like  scheme  has  one  obvious  advantage  over  the 
sodium-neon  approach;  higher  quantum  efficiency.  As  seen  below,  there  are  also 
disadvantages.  It  is  the  purpose  of  this  paper  to  report  a  detailed  analysis  of 
the  A1  XI  -  Mg  IX  resonant  laser,  and  quantify  the  design  requirements  for  a 
demonstration  of  gain. 


II.  ATOMIC  MODELS  AND  CALCULATIONS 

A.  Aluminum  pump  plasma 

The  resonant  pump  line  is  the  Al  XI  ls^2s  ^S2/2“^s^3p  ^Pi/2  transition  at 
0 

48.338  A  (256.49  eV),  nearly  perfectly  coincident  with  the  pumped  Mg  IX  line 

ls^2s^  ^Sq  -  ls^2s4p  at  48.340  A  (256.48  eV).  The  J  =  3/2  component  of  the 

Al  XI  resonance  doublet  is  located  at  48.297  A  (256.71  eV).  This  J  =  3/2  line 

0 

IS  separated  by  43  mA  (8  Doppler  widths  at  an  ion  temperature  of  150  eV)  from 
the  pumped  transition,  but,  if  sufficiently  broadened,  may  participate  in  the 
pumping.  Any  attempt  to  calculate  the  pumping  rate  must  therefore  include 
multi frequency  radiative  transfer  to  resolve  the  doublet  and  obtain 
monochromatic  power. 

The  atomic  model  employed  in  our  calculations  for  aluminum  is  a  subset  of 
the  comprehensive  structure  described  in  Ref.  17.  For  obvious  reasons,  Al  XI 
and  nearby  stages  are  emphasized.  In  the  key  Li-like  Al  XI  species,  the 
following  levels  are  modeled:  (ls^)2s^S,  2p^P,  3s^S,  3p^P,  3d^D,  4s,  4p,  4d, 

4f,  and  a  consolidated  n=5  state.  Note  that  the  J=l/2  and  J=3/2  levels  of  the 
3p^P  level  are  lumped  together;  in  calculating  the  emergent  line  power  profile 
for  this  doublet  the  J-sublevels  are  assumed  to  be  statistically  populated.  The 
Al  plasma  is  assumed  to  be  a  uniform  cylinder  of  specified  density  and 
temperature  in  collisional-radiative  equilibrium  (CRE).  The  populations  are 
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computed  by  solving  the  steady-state  rate  equations  coupling  the  levels, 
including  photon  transport  initially  performed  using  a  multicell  coupling  escape 
probability  technique  described  in  Refs.  18  and  19.  Once  the  populations  are 
computed  by  using  this  efficient  technique,  the  detailed  line  profile  of  the 
pump  doublet  is  obtained  by  a  one-step  multifrequency  recomputatiun  of  the 
radiative  transfer  within  the  two  closely  spaced  lines  assuming  Voigt  profiles. 
Approximately  200  frequency  points  are  employed  in  this  final  phase  of  the 
calculation. 

B.  Pumped  magnesium  plasma 

The  magnesium  atomic  model  has  been  constructed  along  lines  similar  to 
that  of  aluminum,  with  the  obvious  exception  that  the  Be-like  lasing  stage  is 
emphasized.  In  calculating  level  energies,  oscillator  strengths,  and  collision 
rates  for  Mg  IX,  17  different  configurations  are  allowed  to  mix.  The  levels, 
their  statistical  weights,  and  energies  calculated  for  Mg  IX  using  the  atomic 
structure  code  of  Cowan^^  are  presented  in  Table  I.  The  magnesium  plasma  is 
assumed  optically  thin  in  all  cases  considered  below.  The  effects  of  line 
trapping  on  gain  are  neglected  although  this  effect  may  be  non-negligible  at  the 
higher  densities  considered,  for  realistic  plasma  sizes.  Such  an  investigation 
is  beyond  the  scope  of  this  paper  and  is  reserved  for  the  future.  Lines  which 
may  exhibit  gain  with  sufficient  resonant  pumping  of  the  4p  singlet  level  are 
principally  the  4-3  singlet  and  triplet  allowed  transitions.  In  calculating  the 
gain  for  the  various  candidate  lines,  statistical  equilibrium  among  the  singlet 
and  triplet  sublevels  is  assumed.  By  examining  the  results  of  the  gain 
calculations  described  below,  it  is  evident  that  the  most  attractive  laser  line 
by  far  is  the  2s3d  ^D2  -  2s4f  transition  calculated  to  lie  at  228  A.  Other 

candidate  laser  lines  are  predicted  to  have  maximum  gain  more  than  a  factor  of  3 
below  that  of  this  line.  All  further  discussion  of  achievable  gain  therefore 
refers  to  the  2s3d  ^02  -  2s4f  transition. 
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Ill .  RESULTS 


A.  Aluminum  plasma  pump  power 

Fine  aluminum  vires  have  long  served  as  a  source  of  plasma  in  multi- 
megampere  Z  pinch  experiments  (see,  e.g.,  Refs.  21  and  22  and  references 

therein).  However,  by  far  the  most  effort  using  aluminum  has  been  expended  in 
studying  and  analyzing  the  behavior  of  its  K  shell  lines.  Fortunately,  some 
recent  experiments  using  the  3-4  MA  Double  EAGLE  generator  of  Physics 

International  Company^^*^^  provide  considerable  insight  and  guidance  for  the 
optimization  of  L  shell  radiation  which  includes  the  48.3  A  Al  XI  pump  line.  In 

a  series  of  shots  in  which  the  product  of  the  Al  wire  array  mass  (H)  and  the 

square  of  its  initial  radius  (Rq^)  remained  approximately  constant,  the  behavior 
of  L-shell  and  K-shell  radiation  was  observed.  Keeping  M  R^^  constant  (while 
separately  varying  M  and  Rq)  keeps  the  energy  coupling  from  the  generator  to  the 
imploding  load  constant  and  approximately  optimum.  It  was  observed  in  these 
experiments^^’^^  that  the  L  shell  radiation  maximized  at  a  wire  mass  load  of  474 
yg/cm,  and  that  the  diameter  of  the  L  shell  emitting  region  was  approximately 
4.2  mm.  Therefore,  as  a  means  of  assessing  the  achievable  pump  power,  a  series 
of  calculations  was  performed  assuming  a  2-cm  long  cylindrical  aluminum  plasma 
of  diameter  4.2  mm  containing  948  yg.  This  corresponds  to  an  atomic  density  of 
7.6  X  10  ^  cm  The  temperature  was  varied  independently  from  75  to  600  eV  to 

o 

observe  the  effect  on  the  emitted  power  of  the  resonant  pump  line  at  48.338  A 
(256.49  eV)  and  other  lines. 

Results  for  three  temperatures  for  the  Al  XI  2s^Sj/2  “  ^P^^l/3  3/2 

are  shown  in  Fig.  1.  In  the  optically  thin  limit  the  non-pumping  J  =  3/2 

0 

component  at  48.297  A  (256.71  eV),  oscillator  strength  (0.22)  is  twice  as 
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strong  as  the  J  =  1/2  pump  line.  For  the  400  eV  profile  depicted  in  Fig.  1, 
this  is  nearly  the  case.  At  such  a  relatively  high  temperature,  the  Li-like 
stage  is  nearly  burned  through,  and  the  optical  depth  of  the  lines  are  of  order 
unity.  The  two  doublet  components  are  distinct,  with  the  J  =  3/2  line  not  quite 
twice  as  strong  as  the  J  =  1/2.  As  lower  temperatures  are  assumed,  more  Li-like 
A1  results  in  higher  optical  depth.  Since  the  J  =  3/2  component  has  twice  the 
opacity  of  the  J  =  1/2,  the  emitted  intensity  of  the  two  lines  tends  to 
equalize,  as  is  clearly  seen  for  the  250  eV  temperature  case  in  the  center  of 
Fig.  1.  At  still  lower  temperatures  (150  eV),  opacity  broadening,  self 
reversal,  and  the  relatively  lower  optical  depth  between  the  two  components 
produces  the  three-pronged  profile  also  seen  in  Fig.  1.  The  most  important 
implication  of  these  results  for  resonant  photopumping  is  that  monochromatic 
powers  near  100  GW/eV  (10^^  photons  sec“^  Hz“^)  are  achieved  for  a  broad  range 
of  temperatures  (125  -  400  eV).  As  the  Li-like  stage  is  depleted  at  higher 
temperatures,  the  reduced  opacity  of  the  line  keeps  the  pump  power  from  falling 
until  the  line  is  thin  and  there  is  virtually  no  A1  XI  remaining.  In  the  gain 
calculations  described  in  the  next  section,  the  pump  power  (to  4n)  has  been 
fixed  at  100  GV/eV,  a  reasonable  estimate  of  the  capabilities  of  an  optimized  Al 
array  driven  by  several  MA. 

Before  presenting  the  Mg  gain,  one  final  feature  of  these  Al  calculations 

0 

merits  discussion.  Li-like  Al  emits  several  strong  lines  from  33  to  55  A, 
including  the  pump  doublet  discussed  above.  Some  of  these  lines  can  photoionize 
lower  stages  of  Mg,  affecting  the  ionization  balance  of  the  lasant  species  as  a 
function  of  temperature.  Table  II  presents  the  total  power  in  the  more 
important  photoionizing  lines  as  calculated  at  a  temperature  of  125  eV.  Their 
influence  on  Mg  is  discussed  below. 
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B.  Magnesium  gain  without  photoionization 


It  is  assumed  that  a  suitable  Mg  plasma  can  be  placed  in  close  proximity 
to  the  A1  pinch.  In  fact,  separations  of  1-2  cm  in  similar  experiments  have 
been  demonstrated^^.  The  Mg  plasma  may  require  energy  input  other  than  the  x- 
rays  from  the  Al  pinch,  for  instance  a  modest-current  pulser  or  laser  driver. 
In  the  present  work,  we  assume  the  temperatures  and  densities  for  optically  thin 

O 

Mg  plasma  and  calculate  the  gain  in  the  228  A  laser  transition  given  the  Al 
photoresonant  power.  The  pump  rate  from  the  Be-like  ground  state  is 

2 

P(sec“^)  =  f  F  <T  dv  =  F  f  a  dv  =  —  F  f  (1) 

'  ‘'Jvv  vjv  mcv 

where  F^  is  the  monochromatic  pump  flux,  assumed  constant  across  the  pumped  line 
profile,  is  the  photoabsorption  cross-section,  and  f  is  the  absorption 
oscillator  strength  of  the  pumped  transition  (=  0.155).  For  a  2-cm  Al-Mg 
separation,  F^  =  2x10^^  photons  cm~^  sec“^  Hz~^,  for  our  calculated  pump 

strength  of  10^^  photons  sec"^  Hz"^  and  50  dilution.  Using  Eq.  (1)  a  pump 
rate  of  8.2  x  10®  sec“^  is  obtained.  For  a  1  cm  separation,  this  rate  increases 
by  about  a  factor  of  4,  to  3.3  x  10^  sec"^. 

Since  the  Is^  2s^  Be-like  ground  state  forms  the  lower  level  of  the 

pumped  transition,  it  is  important  to  optimize  its  fractional  population  to 

obtain  an  efficient  laser.  Figure  2  displays  the  calculated  Be-like  ground 

state  fraction  as  a  function  of  temperature  and  ion  density.  In  this  case,  the 

only  external  radiation  effect  included  is  the  photoresonant  pumping  at  an 

assumed  separation  of  2  cm.  Inspection  of  Fig.  2  reveals  that  the  maximum 

fractional  population  is  less  than  5%,  near  90  eV  and  10  ions  cm 

0  1  1 

Furthermore,  gain  in  the  228  A  2s3d  ^D2  -  2s4f  -^F^  transition  (not  shown)  does 
not  exceed  0.03  cm“^  anywhere  in  the  densi ty-temperature  plane.  The  key  to  this 
result  is  found  in  the  energy-level  structure  of  Mg  IX,  Table  I.  There  are  5 
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terms  within  61.2  eV  of  the  ground  state  having  a  total  statistical  weight  of 
27,  as  opposed  to  unity  for  the  ls^2s^  ground  level.  Therefore,  it  is 

virtually  impossible  to  obtain  a  large  fractional  ground  state  population. 
Moreover,  the  lower  laser  level  (only  220  eV  from  the  ground)  is  susceptible  to 
collisional  excitation  from  the  ground  state  as  well  as  the  highly  populated 
levels  within  61.2  eV  of  the  ground.  The  fact  that  the  Be-like  ground  state 

maximizes  near  T  =  90  eV  has  the  consequence  that  the  lower  laser  level  is 

easily  populated  via  collisional  excitation.  In  summary,  if  resonant 
photoexcitation  from  the  A1  pinch  is  the  only  external  radiation  considered, 
gain  in  Mg  IX  is  very  low  (<  0.03  cm“^)  due  to  the  small  achievable  population 

of  the  pumped  Be-like  ground  state  as  well  as  strong  collisional  excitation  of 
the  lower  laser  level. 

C.  Gain  calculations  including  photoionization 

Ve  find  that  this  laser  shows  substantially  improved  gain  when 

photoionization  of  Mg  due  to  L-shell  Al  lines  is  considered.  This  allows  the 
Be-like  stage  to  maximize  at  much  lower  temperatures,  reducing  collisional 
excitation  to  the  luwar  laser  level  and  narrowing  the  gain  line  Doppler  width. 
The  ground  state  ionization  potentials  for  Mg  VI,  VII  and  VIII  are, 

respectively,  187,  225,  and  266  eV.  Also  the  ionization  potential  from  the 
critical  2s2p^  first  excited  level  of  B-like  Mg  VIII  is  229  eV.  According  to 
Table  II,  many  of  the  powerful  lines  of  Li-like  Al  are  sufficiently  energetic  to 
ionize  these  stages  toward  the  desired  Mg  IX  ground  state.  The  effect  of  these 
lines  on  the  Mg  ionization  balance  is  revealed  in  the  contour  plots  of  Figs.  3 
and  4,  where  the  Be-like  ground  state  fractional  population  is  plotted  vs.  ion 
density  and  temperature  with  photoionization  included.  Note  that,  for  both  the 
2  cm  separation  (Fig.  3)  and  1  cm  (Fig.  4),  the  Mg  IX  ground  state  fraction 
maximizes  at  much  lower  temperatures,  10-20  eV.  The  effect  on  gain  for  both  the 
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2  cm  case  (Fig  5)  ai.J  1  cm  (Fig.  6)  is  very  positive.  For  the  1  cm  separation 
gains  exceeding  3  cm“^  are  possible  near  ion  densities  of  10  cm“^,  and  for  2 
cm,  the  gain  is  a  more  modest  0.4  cm~^,  peaking  at  somewhat  lower  densities. 
Ion  densities  above  10^®  cm“®  have  not  been  considered  due  to  the  probable 
deleterious  effects  of  line  trapping  for  realistic  plasma  sizes  at  such  higher 
densities. 

The  importance  of  photoionization  of  Mg  by  Al  XI  lines  to  the  overa.^l* 
ionization  balance  is  demonstrated  by  comparing  rates  obtained  from  well-known 
simple  formulations  for  both  photoionization  and  collisional  ionization. 
Consider  the  ionization  from  the  Is  2s2p  levels  of  B-like  Mg  to  the  pumped  Be¬ 
like  ground  state.  The  ionization  potential  of  7.29  eV  pert. .ts  very  efficient 
photoionization  by  the  Al  XI  2p-3d  and  2s -3p  lines  at  237  and  257  eV, 
respectively.  According  to  Table  II,  the  combined  power  of  these  lines  is  185 
GW,  giving  a  flux  at  a  distance  of  i  cm  of  3.7x10^®  ucons  sec“^. 

Multiplying  this  flux  by  the  scre-neu  hydrogenic  Kramers-Gaunt  edge  cross 
section  of  3.6x10“^^  cm  ^  gives  a  photoionization  rate  of  1.3x10®  sec~^.  By 
contrast,  at  45  eV,  the  collisional  ionization  rate  coefficient  obtained  from 
the  Lotz  torn'  .cion  is  4.8  x  10  cm-*  sec"  .  Even  at  the  highest  electron 
densities  consiJered  here,  10^^  cm”®,  the  corresponding  rate  of  4.8x10^  sec”^  is 
nearly  a  factor  of  3  below  that  of  photoionization.  At  lower  temperatures  and 
den.'-^ties,  the  disparity  is  much  greater. 

The  effect  of  ionization  and  collisional  excitation  on  gain  is  illustrated 
by  close  study  of  Fig  6  at  an  ion  density  of  3.6x10  cm  .  At  a  temperature  of 
12  eV,  gain  of  2  cm"^  is  predicted,  whereas  at  65  eV,  only  0.045  cm”^  is 
expected.  The  fractional  population  of  the  upper  laser  level  is  given  within  a 
factor  of  2  under  these  conditions  by  a  simple  kinetics  model 

f(2s  4f  ^F3)  =  f(g)  •  P  •  (7/16)  /  A  (2) 


where  f(g)  is  the  pumped  Be-like  ground  state  fraction,  P  is  the  photoresonant 


pump  rate  (3.3  x  10  ^  sec“^),  7/16  is  the  statistical  weight  of  the  4f  level 
divided  by  that  of  all  of  the  2s4f  singlet  levels,  and  A  is  the  weighted  inverse 
lifetime  of  the  n=4  singlets  (1.6  x  10^^  sec“^).  Note  from  Fig. 4  that  f(g)  is 
reduced  from  6.5  x  10“^  to  6.2  x  10“^  as  the  temperature  increases  from  12  to  65 
eV.  According  to  Eq.  (2),  this  would  reduce  the  gain  from  2  cm“^  to  0.19  cm“^. 
The  increase  in  line  Doppler  width  accounts  for  a  further  reduction  of  a  factor 
of  2.3.  The  final  factor  of  2  reduction  required  arises  from  the  increase  in 
lower  laser  level  population.  At  12  eV,  the  4f  population  exceeds  that  of 
3d  ^D2  by  a  factor  of  5.6;  at  65  eV  this  ratio  is  reduced  to  2.2  (as  given  by 
the  detailed  calculations).  This  is  not  surprising  as  the  collisional 
excitation  rate  from  the  Be-like  ground  state  increases  5  orders  of  magnitude 
between  12  and  65  eV.  The  above  simple  considerations  provide  a  useful  guide 
and  scoping  tool  for  the  kinetics  as  revealed  by  the  detailed  calculations. 

IV.  CONCLUDING  DISCUSSION 


Our  detailed  calculations  of  photoresonant  pump  power  in  A1  XI  as  well  as 
gain  in  Mg  IX  reveal  that  photoionization  by  nonresonant  Al  lines  is  required 
along  with  the  resonant  photopumping  to  obtain  useful  gain.  The  photoionization 
allows  the  Mg  IX  to  exist  at  relatively  low  electron  temperatures  greatly 
reducing  deleterious  collisional  excitation  to  the  lower  laser  levels.  It  is 
instructive  to  compare  this  scheme  to  the  well-known  photoresonant  Na  X  -  Ne  IX 
system. ^’®“^®Such  a  comparison  is  clarified  and  facilitated  by  the  fact  that  the 
gain  lines  have  nearly  the  same  wavelength  and  oscillator  strength-indeed ,  they 
are  both  3d  ^D2-4f  ^Fj  transitions.  However,  the  quantum  efficiency  of  the  Al- 
Mg  system  exceeds  that  of  the  Na-Ne  scheme  by  a  factor  of  4.4.  That  is,  one 
pumping  of  the  Ne  IX  upper  laser  level  requires  a  photon  4.4  times  more 
energetic  than  that  required  to  pump  Mg  IX.  Moreover,  this  clear  advantage  for 
Al-Mg  is  offset  by  the  difficulty  of  obtaining  a  large  fraction  of  Mg  IX  ions  in 
the  pumped  ls^2s^  ground  state,  due  to  the  proximity  of  many  other  levels  of 
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much  higher  statistical  weight.  By  contrast  the  Ne  IX  pump  ground  level  (Is 
^Sq)  is  isolated  from  all  overlying  levels  by  nearly  a  kilovolt  and  at  least  50% 
of  the  neon  ions  can  be  placed  in  this  state  in  a  properly  prepared  plasma.  The 
end  result,  even  when  many  other  atomic  physics  details  are  included,  is  that 
the  two  systems  are  about  equally  attractive. 
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TABLE  I:  Levels,  their  statistical  weights,  and  energies  in  the  Be-like 

magnesium  atomic  model. 


Level 

Statistical  Weight 

Energy  (eV) 

Ishs^ 

1 

0. 

ls^2s2p  3p 

9 

17.8 

Is22s2p 

3 

32.4 

Is22p2  3p 

9 

45.6 

Is22p2  Id 

5 

50.2 

ls^2p2 

1 

61.2 

1s22s31  ^SPD 

27 

199.3 

1s22s31  ^SPD 

9 

201.1 

Is22p31  ^SPDF 

81 

220.4 

ls^2p31  ^SPDF 

27 

222.8 

ls^2s41  ^SPDF 

48 

257.8 

15^2841  ^SPDF 

16 

258.4 

15^2551 

100 

283.7 

ls^2s61 

144 

297.4 
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TABLE  II;  Total  power  output  of  various  photoionizing  transitions  in  Li-like 
aluminum.  Length  of  Z  pinch  is  2  cm,  diameter  4.2  mm,  ion  density  7.6xl0^^cm“^, 
and  temperature  125  eV.  Doublets  and  other  closely  spaced  lines  are  lumped 
together. 


Transition  (Is^  core)  Energy  (eV) 

2p-3s  228 

2p-3d  237 

2s-3p  257 

2p-4s  314 

2p-4d  317 

2s-4p  338 

2p-51  354 

2s-51  376 


Power  (GW) 

28 

98 

87 

20 

143 

76 

19 

14 


28 


HOWER  (GW/eV) 


Fig.  1. 


400eV 


256  257 


256  257 


ENERGY  (eV) 


o  2 

Power  output  line  profile  for  the  Al  XI  2s  ^^\/2  ~  ^1/2, 3/2  P'^'^P 

line  doublet  for  an  assumed  2  cm  long  Z  pinch  with  ion  density  7.6  x 
10^^  cm“^.  Temperatures  for  the  3  cases  are  indicated.  The  pumped  Mg 
transition  lies  at  256.48  eV. 


10'^  10^® 

ION  DENSITY  (cm -2) 

Be-like  ground  state  fraction  for  an  optically  thin  Mg  plasma  vs.  ion 
density  and  temperature.  Distance  of  Mg  plasma  from  Al  Z  pinch  is 
2  cm.  Resonant  pumping  of  Mg  IX  2s4p  level  is  included,  but  not 
photoionization  of  Mg  by  the  Al  radiation. 


Temperature  (eV) 


Ion  Density(cm"^) 


Fig.  3.  As  in  Fig.  2,  except  that  photoionization  of  the  various  Mg  stages 
the  A1  lines  listed  in  Table  II  is  included. 
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Fig.  5.  Gain  coefficient  for  the  Mg  IX  2s3d  ^02  -  2s4f  ^F3  transition  near 
228  A  is  shown  as  a  function  of  ion  density  and  temperature  for  an 
Al-Mg  separation  of  2  cm,  with  both  photoionization  and  resonant 
photopumping  included. 
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Design  Considerations  for  Z-Pinch  Driven  Photoresonant 
X-Ray  Lasing  in  Neonlike  Krypton 

J.  V.  Thornhill,  J.  P.  Apruzese,  J.  Davis,  and  R.  U.  Clark 
Naval  Research  Laboratory 
Radiation  Hydrodynamics  Branch 
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Washington,  D.C.  20375-5000 

Abstract 

A  neonlike  x-ray  laser  photoresonant  pumping  scheme  is  explored.  An 
attractive  design  is  a  coaxial  z-pinch  consisting  of  an  inner  krypton  lasant 
plasma  surrounded  by  a  carbon  shell  that  itself  is  surrounded  by  a  stagnated 
krypton  plasma.  The  photoresonant  radiation  emitted  from  the  outer  plasma 
passes  inwardly  through  the  carbon  shell  and  photoexcites  lasant  electrons 
to  the  3s,  3d,  and  4d  levels.  It  is  calculated  that  monochromatic  pump 
powers  of  30  GW/cm-eV  can  be  achieved  for  the  4d  neonlike  resonance  line 
while  powers  of  100  -  200  GW/cm-eV  are  attainable  for  the  3s  and  3d 

resonance  lines.  The  gain  in  several  neonlike  3s-3p  transitions  is 

calculated  as  a  function  of  temperature  and  density  of  the  lasant  plasma. 
Reasonable  gain  in  only  the  3p-3s  (J=0-1)  transition  is  obtained  for  high 

density,  high  temperature  lasant  conditions.  This  gain,  which  is  driven  by 
monopole  excitation  from  the  ground  state  to  the  3p  J=0  level,  is  nearly 
independent  of  the  presence  of  pump  radiation.  For  low  density,  low 

temperature  lasants,  optimal  gain  is  obtained  in  the  3p-3s  (J=2-l) 

transitions.  Under  these  conditions,  the  pump  radiation  is  necessary  to 
photoionize  and  photoexcite  the  lasant  plasma  to  the  neonlike  stage  while 
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radiative  decay  from  the  resonantly  pumped  4d  level  to  the  3p  (J=2)  states 
provides  the  major  upper  lasing  level  pump  mechanism. 

I.  INTRODUCTION 

One  of  the  major  limitations  of  photoresonant  laser  schemes  is  the 
wavelength  requirement  that  the  pumping  and  pumped  transitions  lie  within  a 
linewidth  of  one  another.  There  are  a  number  of  candidate  resonant  pairs 
that  meet  this  requirement.^”^  However,  an  assured  way  of  achieving  this 
goal  is  to  use  ions  of  the  same  kind  for  both  the  pump  source  and  lasant. 
Recent  work^  ^  has  shown  that  the  potential  of  photoresonant  pumping  by  like 
ions  to  produce  gain  in  the  x-ray  region  compares  favorably  to  that  of  other 
photoresonant  schemes  as  well  as  to  collisional  excitation  or  recombination 
schemes  driven  by  high  powered  lasers  or  z-pinches. 

In  this  investigation  we  build  upon  the  analytical  work  of  Chichkov  and 
Fill^”^  and  discuss  the  results  of  a  numerical  model  of  their  proposed 
coaxial  z-pinch  configuration  laser  in  which  the  inner  core  lasant  plasma  is 
being  pumped  by  a  stagnating  outer  plasma,  consisting  of  the  same  element  as 
the  core.  In  particular,  a  neonlike  krypton  photoresonant  scheme  is 
analyzed.  In  the  absence  of  photoresonant  pumping  the  neonlike  3p  (J=0) 
level  is  populated  primarily  by  direct  monopole  excitation  of  the  2p^  ground 
state,  see  Fig.  1,  whereas  the  other  3p  levels  are  populated  by 
recombination  from  the  ground  state  of  fluorinelike  krypton  as  well  as  d  - 
excitation  from  the  upper  levels-  The  presence  of  pump  radiation  produces 
two  main  effects  which  alter  the  populations  of  the  3p  states.  First, 
enhanced  photoexcitation  and  photoionization  ionize  low  density  plasmas  to 
the  neonlike  state  at  much  lover  temperatures  than  are  required  without 
photo-pumping.  Also,  the  photoresonant  pumping  of  the  4d  level  with 
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subsequent  radiative  decay  to  the  3p  states  provides  for  a  strong  upper 
(J*2)  lasing  level  pump  mechanism.  As  mentioned  by  Chichkov  and  Fill,^”^  the 
advantage  of  this  scheme  is  that  it  operates  at  electron  densities  too  low 
to  cause  collisional  de-excitation  of  the  laser  transition  or  refraction  of 
the  amplified  beam.  In  addition,  the  quasi-metastability  of  the  upper  lasing 
level  provides  inversion  without  filtering  of  the  3s  resonance  radiation, 
which  populates  the  lower  lasing  level. 

The  z-pinch  model  developed  to  study  this  scheme  has  been  simplified  by 
assuming  that  the  pump  plasma  has  already  stagnated  against  a  cold  carbon 
shell  that  surrounds  and  isolates  the  core  lasant  plasma  from  hydrodynamic 
effects.  The  radiation  from  the  outer  plasma  passes  inwardly  through  the 
massive  carbon  shell  and  photoexcites  and  photoionizes  electrons  of  the  core 
plasma.  It  is  assumed  that  the  core  plasma  can  be  independently  heated 
either  Ohmically,  by  division  of  the  current  driving  the  z-pinch,  or  else  by 
a  laser  combined  with  radiation  from  the  pump  plasma. 

The  details  of  the  z-pinch  configuration  and  the  multi-material 
radiation  hydrodynamics  models  used  in  the  simulation  are  discussed  in 
Section  II.  Section  III-A  describes  the  calculation  of  the  pump  powers  and 
section  III-B  presents  results  for  the  temperature  and  density  dependence  of 

the  effective  charge  and  gain  in  the  neonlike  3p-3s  (J=0-1)  transition  at 

0  0 
172  A  and  the  (J=2-l)  transition  at  186  A  (Fig.  1).  The  results  for  gain 

0 

achieved  in  the  (J=2-l)  184  A  transition  are  not  discussed,  but  they  are 

O 

very  similar  to  the  results  for  the  (J=2-l)  186  A  transition.  A  contrast  of 
lasant  characteristics  with  and  without  the  influence  of  the  pump  plasma  as 
well  as  a  discussion  of  the  role  that  opacity  plays  in  limiting  the 
achievable  gain  are  also  provided  in  Section  III-B.  A  summarizing 
discussion  follows  in  Section  IV. 
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II.  MODEL  DESCRIPTION 


This  model  describes  a  static  situation  that  is  assumed  to  be 
representative  of  the  conditions  realized  during  the  stagnation  phase  of  a 
2-pinch  experiment.  In  our  model  the  cylindrical  plasma  is  divided  into 
three  regions  (Fig.  2).  The  first  region  contains  the  lasant  krypton  plasma 
and  it  extends  from  the  origin  to  0.1  cm.  It  has  uniform  density  and 
temperature  profiles.  Since  the  lasant  plasma  is  independently  heated,  its 
temperature  is  freely  specified.  Region  2  is  a  thin  column  extending  from 
0.1  cm  to  0.10024  cm,  it  consists  of  400  ug/cm  of  solid  state  density  carbon 
at  5  eV.  If  this  were  a  dynamic  Z-pinch  simulation,  the  purpose  of  the 
carbon  in  hydrodynamically  isolating  the  core  plasma  could  be  tested. 
However,  since  this  calculation  is  static,  it  is  assumed  that  the  carbon  has 
done  its  job  and  that  the  outer  krypton  pump  plasma  has  stagnated  onto  the 
carbon  shell.  At  this  stage  the  main  effect  of  the  carbon  is  to  somewhat 
attenuate  via  inner-shell  absorption  the  radiation  field  emanating  from  the 
stagnated  outer  plasma.  The  third  region,  which  extends  from  0.10024  cm  to 
0.15  cm,  contains  a  krypton  pump  plasma  that  has  uniform  temperature  and 
density  profiles  characteristic  of  a  stagnated  plasma. 

The  analysis  proceeds  in  two  stages.  First,  a  multi-level  non-LTE 
ionization  and  radiation  transport  calculation  that  self-consistently 
couples  all  of  the  regions  of  the  z-pinch  configuration  is  performed.  This 
is  an  extensive  model  that  is  briefly  described  in  Appendix  I.  Because  it  is 
too  expensive  computationally,  there  is  insufficient  neonlike  krypton  level 
structure  represented  in  this  model  to  do  a  credible  laser  gain  calculation. 
Instead,  the  gain  calculation  is  performed  by  passing  the  following  lasant 
plasma  data  along  to  a  more  detailed  multi-level  gain  model  that  completes 
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the  second  stage  of  the  calculation,  this  data  includes;  the  electron  and 
ion  densities,  ground  state  populations  of  neonlike  and  fluorinelike 
krypton,  and  the  pump  powers  of  the  3s  (3/2, 1/2)^,  3s  (1/2, 1/2)^,  3d 
(3/2,5/2)j^,  3d  (l/2,3/2)j^,  and  4d  neonlike  resonance  lines  I(j2^,j2)j 
notation]  emanating  from  the  outer  plasma,  see  Fig.  1.  This  gain  model  is 
similar  to  the  neonlike  selenium  model  described  in  Ref.  7,  but  scaled  to 
neonlike  krypton.  A  static  local  approximation  is  used  to  transport 
radiation.  The  pump  powers,  attenuated  by  e  v  ,  where  is  the  frequency- 
dependent  optical  depth  and  x  is  the  distance  into  the  core  plasma,  are 
included  as  source  terms  in  the  gain  model  for  populating  the  3s,  3d,  and  4d 
levels. 

Ill .  RESULTS 


A.  Krypton  plasma  pump  power 

At  stagnation  the  krypton  pump  plasma  is  specified  to  have  a  mass 

19  -3 

loading  of  273  ug/cm,  an  ion  density  of  5.0  x  10  cm  ,  and  temperature  of 
500  eV.  Good  pumping  rates  over  a  wide  range  of  temperatures  and  densities 
were  found.  The  chosen  conditions  represent  typical  values  that  are  believed 
to  be  attainable  (similar  values  have  been  obtained  for  nickel  wire 

g 

experiments  performed  at  Physics  International  Corp.)  using  existing  pulse 
power  machines,  i.e,  Sandia  National  Laboratories'  Saturn  machine,  Physics 
International's  Double  Eagle  Machine,  or  Maxwell  Laboratories'  Blackjack  5 
device. 

The  powers  of  the  3s,  3d,  and  4d  resonance  lines  of  the  pump  plasma  are 
shown  in  the  theoretical  spectrum  displayed  in  Fig.  3.  For  illustration 
purposes  an  enlargement  of  the  3s  (3/2, 1/2)^  resonance  line  is  shown  in  Fig. 
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4.  According  to  this  figure,  monochromatic  powers  near  170  GW/eV-cm  are 
achievable  for  the  3s  (3/2, 1/2)^^  resonance  transition.  Similar  to  the 
calculation  of  Ref.  9,  the  pump  rate  for  the  3s  (3/2, 1/2)^  level  (for 
example)  from  the  neonlike  ground  state  is  calculated  as 

P(sec'^)  =  J  dv  =  J  dv  =  (ne^/mc)  F^  f  (1) 

where  F^  is  the  monochromatic  pump  flux,  assumed  constant  across  the  pumped 

line  profile,  is  the  photoabsorption  cross-section,  and  f  is  the 

absorption  oscillator  strength  of  the  pumped  transition  (=0.136).  Assuming 

negligible  attenuation  of  the  pump  radiation  by  the  carbon  shell  and  that 

the  core  plasma  sees  the  full  2n  steradians  of  radiation  field  from  the 
12  2 

pump,  F^  =  2.7x10  photons/(cm  -sec-Hz).  Using  Eq.  (1)  a  pump  rate  of  9.7  x 
9  -1 

10  sec  is  obtained.  The  pump  rates,  energy  levels,  and  oscillator 

strengths  for  the  other  3s,  3d,  and  4d  levels  are  presented  in  Table  I.  The 

pump  rate  to  the  4d  state,  which  ultimately  populates  the  upper  3p  (J=2) 

9  -1 

lasing  levels  by  radiative  decay  is  8.29  x  10  sec  .  Whether  this  factor  is 
sufficient  to  produce  or  enhance  gain  is  examined  in  the  next  subsection. 

B.  Calculation  of  gain  in  krypton 

The  effective  charge  as  a  function  of  temperature  and  density  of  the 
lasant  plasma  for  the  case  in  which  the  radiation  field  emitted  from  the 
pump  plasma  is  turned  off  is  shown  in  Fig.  5.  These  results  are  taken  from 
the  innermost  zone  of  the  lasant  plasma,  but,  they  also  represent  typical 
values  throughout  the  core.  This  figure  illustrates  that  increasing  the  core 
ion  density  allows  the  neonlike  krypton  state  (z  =  25)  to  be  attained  at 
progressively  lower  temperatures.  The  principal  cause  of  this  behavior  is 
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ladder  ionization.  As  density  increases,  the  increased  excited  state 

populations  are  ionized  to  the  next  highest  ground  stage,  requiring  a  lover 

ionization  potential  than  ground-to-ground  ionization. 

Vhen  the  pump  plasma  is  present,  photoionization  and  photoexcitation 

produce  a  dramatic  rise  in  the  ionization  state  of  the  low  density  lasant 

plasmas.  This  is  clearly  seen  in  Fig.  6,  where  the  effectivt  charge  as  a 

function  of  density  and  temperature  is  shown.  A  comparison  of  Fig.  5  and 

Fig.  6  reveals  that  the  low  density,  =10  cm  ,  plasma  increases  from 

an  ionization  stage  of  z  =  22  (without  pump)  to  z  =  26.5  (with  pump)  at  a 

temperature  of  only  250  eV.  In  fact.  Fig.  6  demonstrates  that  the 
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simulations  for  the  10  cm  plasma  should  be  performed  at  lower 

temperatures  or  lower  pump  powers  because  the  neonlike  krypton  stage  is 
partially  burned  through  at  250  eV.  Low  temperature  and  density  results  are 

discussed  below.  The  contrasting  behavior  between  the  state  of  the  lasant 

plasma  with  and  without  a  pump  field  is  an  excellent  example  of  the  trade¬ 
off  between  collisional  versus  radiative  effects,  i.e.,  the  state  of  low 
density  (less  collisional)  plasmas  are  significantly  affected  by  the  pump 
field,  whereas,  higher  density  (collisional)  plasmas  are  only  minimally 
altered. 

When  there  is  no  pump  plasma,  significant  gain  is  only  achieved  for 
19  -3 

the  high  density,  10  cm  ,  plasmas  principally  because  only  they  attain 

sufficient  population  in  the  excited  neonlike  states,  see  Fig.  5.  The  gain 

0 

in  the  3p  -  3s  (J=0-1)  transition  at  172  A  is  plotted  in  Fig.  7  for  the 

inner  zone  of  the  lasant  plasma  as  a  function  of  temperature.  This  gain  is 

nearly  uniform  from  the  center  to  within  100  I'm  of  the  edge  of  the  lasant 
plasma.  At  the  very  edge  of  the  lasant  plasma  reduced  opacity  of  the  3s 
(\/2,l/2)^  resonance  line  does  lead  to  factors  of  2  -  3  enhancement  of  gain 

over  that  displayed  in  Fig.  7.  Gain  in  excess  of  unity  was  achieved  in  the 
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(J«2-l)  transition  only  at  the  very  edge  of  the  plasma.  Again,  this 
phenomenon  is  the  result  of  reduced  trapping  of  the  3s  (1/2, 1/2 )j^  resonance 
line.  This  gain  is  not  shown,  but  it  has  a  value  of  2.25  cm"^  at  700  eV  and 
linearly  extrapolates  to  zero  at  a  temperature  of  450  eV. 

These  static  model  results  are  in  reasonable  accord  with  predictions 
from  the  Livermore  model^^  that  was  used  to  simulate  neonlike  selenium  laser 
experiments.  In  particular,  both  models  show  that  monopole  excitation  of  the 
upper  3p  J=0  level  exceeds  the  indirect  excitation  of  the  3p  J=2  levels  such 

O 

that  gain  in  the  neonlike  selenium  183  A  (J=0-1)  transition  is  substantially 

o 

larger  than  in  the  209  A  (J=2-l)  transition.  These  model  results  are  in 

opposition  to  the  experimentally  observed  amplifications  of  5  cm~^  for  the 

(J=2-l)  transition  and  negligible  gain  in  the  (J=0-1)  transition. This 

'mystery  of  the  missing  (J=0-1)  line'^^  has  led  to  some  controversy  as  to 

whether  the  Livermore  experiments  achieved  the  predicted  plasma  conditions 

necessary  to  produce  high  gain  in  the  (J=0-1)  transition  and  whether  the 

atomic  models  are  flawed.^  This  controversy  has  somewhat  abated  due  to  new 

experiments  which  showed  2-3  cm“^  gain  in  the  (J=0-1)  transition  of 
12 

selenium.  In  addition,  the  Livermore  model  has  been  updated  by  including 

detailed  dielectronic  recombination  to  the  excited  states  of  neonlike 
13-17 

selenium  (with  the  exception  of  an  effective  dielectronic  recombination 

rate  to  the  ground  state,  these  individual  processes  are  not  included  in  our 

model).  This  addition  has  increased  the  predicted  strengths  of  the  (J=2-l) 

lines  relative  to  the  (J=0-1)  line.  However,  the  Livermore  model  still 

12 

predicts  larger  gain  in  the  (J=0-1)  transition.  Another  CRE  selenium  model 

18 

that  includes  dielectronic  recombination  to  the  excited  states  as  well  as 

excitation  and  ionization  coefficients  derived  from  the  same  data  base  as 

the  Livermore  model  has  confirmed  the  prediction  that  the  (J=0-1)/(J=2-1) 

0 

gain  ratio  is  larger  than  unity  for  the  183  and  209  A  transitions  but  it 
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predicts  the  ratio  to  be  less  than  unity  for  the  (J=0-1)  169  A  and  (J=2-l) 

O 

206  A  transitions.  This  atomic  model,  the  Livermore  model,  and  the  one  used 
in  our  work  all  have  common  features  and  well  as  some  that  are  still 
controversial.  These  unresolved  experimental,  atomic  data,  and  calculational 
issues  do  not  affect  photon  pump  rates  and  therefore  do  not  negate  the 
utility  of  evaluating  the  effects  of  photoresonant  pumping  on  a  neonlike 
gain  calculation. 

The  gain  results  are  now  discussed  for  the  case  in  which  the  pump  is  on 
but  the  explicit  pump  rates  to  the  3s,  3d,  and  4d  levels,  listed  in  Table  I, 
are  unattenuated  as  the  pump  radiation  travels  into  the  core.  The  enhanced 
pumping  of  the  3p  J=0  state  by  direct  radiative  decay  of  the  4d  level  does 
lead  to  slightly  larger  (J=0-1)  gain  coefficients  than  displayed  in  Fig.  7 
for  the  non-photo-pumped  lasant  plasma.  The  additional  ionization  and 
enhanced  upper  laser  level  pumping  associated  with  the  pump  plasma  also 
produce  unity  (J=0-1)  gain  coefficients  for  temperatures  in  excess  of  350  eV 
in  moderate  density,  10  cm"  ,  plasmas.  However,  the  most  dramatic  increase 

O 

in  gain  occurs  for  the  (J=2-l)  186  A  transition.  This  is  illustrated  in  Fig. 

8  where  excellent  gain,  >  2  cm"^,  is  achieved  for  temperatures  in  excess  of 

300  eV.  There  are  several  phenomena  that  play  a  role  in  producing  gain 

under  these  plasma  conditions.  First,  since  opacity  scales  with  ion  density, 

the  trapping  of  the  3s  to  ground  state  resonance  lines  is  substantially 
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lower  for  a  plasma  with  an  ion  density  of  10  cm  as  opposed  to  10  cm 
Secondly,  pumping  excites  the  plasma  to  the  neonlike  stage  at  lover 
temperatures  than  without  the  pump,  see  Figs.  5  and  6;  this  reduces  the  line 
Doppler  width  and  also  promotes  recombination  from  the  fluorinelike  stage. 
Third,  radiative  decay  of  the  resonantly  pumped  4d  level  is  the  dominant 
mechanism  for  populating  the  upper  lasing  3p  J=2  level. ^  ^  The  presence  of 
resonant  photo-pumping  to  the  4d  level  effects  a  three  order  of  magnitude 
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enhancement  of  the  creation  rate  to  the  3p  J=2  level.  For  the  10  cm~ 

O 

plasma,  (J=2-l)  186  A  gain  is  limited  by  low  neonlike  ground  state 
populations  due  to  burn  through,  i.e.  the  simulations  were  performed  at  too 
high  lasant  temperatures. 

Gain  results  from  the  calculations  that  include  the  attenuation  of  the 
pump  rates  due  to  resonance  absorption  in  the  core  plasma  are  now  examined. 
These  calculations  are  the  most  representative  of  an  actual  z-pinch 

configuration  soft  x-ray  laser.  Since  the  opacity  of  the  4d  resonance  line 

■  -1  -1  18 

IS  greater  than  60  cm  (200  cm  )  in  a  plasma  having  an  ion  density  of  10 

-3  19  -3 

cm  (10  cm  )  and  a  temperature  in  excess  of  400  eV  ,  it  is  not 
surprising  that  at  sufficiently  high  ion  density  the  upper  lasing  level  pump 
rates  are  attenuated  a  short  distance  into  the  core  plasma.  The  attenuation 
is  sufficient  that  the  gain  results  for  the  (J=0-1)  transition  are  virtually 
identical  to  the  case  without  any  pump  plasma  (Fig.  7),  except  at  the  very 
edge  of  the  plasma  where  the  4d  resonance  line  is  unattenuated.  This  point 

e 

is  also  clearly  demonstrated  by  comparing  the  186  A  (J=2-l)  gain  results  at 

the  center  of  the  outer  zone  (Fig.  9),  which  is  only  67  gm  from  the  edge  of 

the  lasant  plasma,  with  the  unattenuated  results  (Fig.  8).  The  gain  at  67  um 

is  reduced  by  attenuation  to  half  the  value  found  at  the  very  edge  of  the 
18  -3 

10  cm~  lasant  plasma.  Attenuation  is  so  complete  at  the  very  center  of 

18  -3  ° 

the  10  cm”  lasant  plasma  that  there  is  no  gain  in  the  186  A  transition. 
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For  the  10  cm  plasma,  the  ion  density  is  low  enough  that  attenuation  of 
the  pump  only  produces  a  40  percent  decrease  in  gain  from  the  outer  compared 
to  the  inner  zone.  This  is  encouraging  because  by  choosing  a  low  lasant 
density  and  temperature  (Fig.  10),  gain  coefficients  in  excess  of  unity  over 
a  large,  1  mm,  radial  expanse  are  obtainable.  At  these  temperatures  the 
radiation  from  the  pump  source  alone  may  be  sufficient  to  heat  the  lasant 
plasma.  The  major  drawback  to  using  a  low  temperature  lasant  plasma  is  that. 
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as  revealed  by  a  time  dependent  ionization  dynamics  calculation,  the  pump 

plasma  needs  to  be  present  for  prohibitively  long,  100  ns,  time  periods  in 

17  -3 

order  to  photoionize  a  plasma  with  ion  density  of  10  cm  and  temperature 
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of  50  eV  to  the  neonlike  state.  By  contrast,  a  10  cm  plasma  at  300  eV 
requires  only  about  20  ns  to  reach  the  neonlike  stage.  These  photoionization 
times  imply  that  the  low  density  gains  shown  in  Figs.  8  and  9  are 
conservative  because  there  will  be  a  substantial  interval  of  time  during 
photoionization  when  conditions  will  be  more  favorable  for  gain  than 
indicated  by  the  steady  state  effective  charge  of  ~  26.5  shown  in  Fig.  6. 


IV.  SUMMARY  AND  DISCUSSION 

This  work  was  intended  to  be  a  preliminary  assessment,  not  a  full 

optimization  determination,  for  investigating  the  potential  of  realizing 

lasing  via  neonlike  krypton  photoresonant  excitation  in  which  ions  of  the 

same  kind  serve  as  both  the  lasant  and  pump  source.  The  scheme  is  for  a 

coaxial  z-pinch  configuration  in  which  the  inner  core  lasant  is  pumped  by  a 

stagnated  outer  shell  plasma.  Gain  coefficients  as  a  function  of  density  and 

temperature  of  the  lasant  plasma  as  well  as  whether  or  not  the  pump  plasma 

0 

is  present  are  calculated  for  the  3p-3s  {J=0-1)  transition  at  172  A  and  the 

O 

3p-3s  (J=2-l)  transition  at  186  A. 

It  is  calculated  that  monochromatic  powers  of  30  GU/eV-cm  for  4d 
resonance  radiation  are  achievable  by  stagnating  a  pump  plasma  onto  a 
massive  carbon  shell  that  hydrodynamically  isolates  the  lasant  plasma  but  is 
transparent  to  the  resonant  radiation.  This  radiation  pumps  the  4d  states  of 
the  lasant  plasma.  The  subsequent  radiative  decay  from  the  4d  to  the  3p 
level  provides  for  the  major  upper  lasing  level  pump  mechanism  for  the  J=2 
states.  Although  the  J=0  level  population  can  be  enhanced  by  radiative  decay 
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of  the  pumped  4d  level,  under  most  (J=0-1)  gain  conducive  plasma  conditions, 
the  most  important  J=0  level  populating  mechanism  is  direct  monopole 
excitation  from  the  ground  state. 

Significant  gain  coefficients  for  the  3p-3s  (J=0-1)  transition  were 

19  -3 

found  only  for  high  ion  density,  10  cm  ,  and  high  electron  temperature  (> 
500  eV)  lasant  plasmas.  Under  these  conditions,  gains  of  the  order  2  cm~^ 
were  present  throughout  the  interior  of  the  lasant  plasma;  at  the  edge  of 
the  plasma  the  gain  was  a  factor  of  2-3  larger  because  of  reduced  opacity  in 
the  lower  lasing  level  3s  resonance  transition.  Because  of  the  dominance  of 
collisional  processes  in  these  high  density  plasmas,  the  gain  was  nearly 
independent  of  the  presence  of  pump  radiation. 

O 

For  high  density  lasant  plasmas,  opacity  effects  limit  gain  in  the  186  A 
3p-3s  (J=2-l)  transition  to  only  the  outer  edge  of  the  plasma.  For  moderate 
and  low  density  plasmas  the  pump  field  is  essential  for  photoionizing  the 
plasma  to  the  neonlike  state  as  well  as  resonantly  photopumping  the  4d  level 
populations.  Radiative  decay  of  the  4d  to  the  3p  J=2  state  is  the  principal 
creation  mechanism  for  the  J=2  upper  lasing  level.  Attenuation  of  the 
photoresonant  4d  pump  radiation  limits  significant  gain  to  only  the  outer 
100  pm  of  a  lasant  plasma  with  an  ion  density  of  10  cm  .  The  gain  at  this 
density  for  temperatures  of  300  -  700  eV  is  between  1-  2  cm  except  at  the 
very  edge  of  the  plasma  where  it  is  factors  of  2-3  larger  because  of 

negligible  attenuation  of  the  pump  radiation  as  well  as  reduced  opacity.  At 

17  -3  -1 

lower  densities,  10  cm  ,  gains  of  the  order  1  -  2  cm  are  achieved  over 

a  broad  radial  expanse  at  very  low  lasant  temperatures,  of  the  order  50  eV. 

The  disadvantage  of  the  low  density  lasant  plasma  is  the  pump  plasma  has  to 

provide  photoionizing/pumping  radiation  for  the  long  time  span  needed  to 

ionize  the  lasant  plasma  to  the  neonlike  state.  The  advantages  of  the  low 
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density  lasant  plasma  are  the  radiation  from  the  pump  source  alone  may  be 
sufficient  to  heat  the  lasant  plasma,  and  refraction  effects  are  minimized. 
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Appendix  I 


The  ionization  dynamics  is  calculated  using  the  multi-material 
collisional-radiative-equilibrium  (CRE)  radiation  and  atomic  physics  models 
developed  at  In  addition  to  all  of  the  ground  states  of  carbon 
and  krypton,  there  are  136  excited  levels  of  krypton  represented  by  the 
model.  There  is  not  any  excited  state  structure  of  carbon  included  in  the 
model  since  the  neonlike  resonance  lines  of  krypton  lie  above  the  carbon  K 
edge,  thus,  neglecting  excited  state  structure  should  produce  an 
overestimate  of  the  inner-shell  absorption  and  result  in  a  conservative  gain 
calculation. 

Since  optical  pumping  by  such  mechanisms  as  photoionization  and 
photoexcitation  allow  plasmas  separated  in  space  to  affect  each  other's 
level  populations,  these  processes  must  be  self-consistently  modeled  by 
transport  of  the  radiation  field.  In  this  calculation,  a  multi-frequency 
transport  model  is  employed  to  calculate  the  level  populations  of  the  core 
plasma  as  they  are  redistributed  by  the  pump  plasma's  radiation  field.  This 
model  transports  1400  frequency  bins  that  span  the  total  frequency  space 
over  which  the  bound-bound,  bound-free,  and  free-free  radiation  is  emitted 
from  the  krypton  and  carbon  plasmas.  The  model  includes  58  spectral  lines 
mostly  from  the  fluorinelike,  neonlike,  and  sodiumlike  krypton  ionization 
stages.  A  Voigt  line  profile  is  assumed,  and  there  are  12  frequency  bins 
resolving  each  line.  In  addition  there  are  2  frequency  bins  designated  for 
each  bound-free  edge.  Transporting  the  radiation  required  zoning  the  plasma 
regions;  in  our  model  we  used  4  zones  to  model  the  core,  3  zones  to  model 
the  carbon  shell,  and  4  zones  to  model  the  pump  plasma.  All  the  zones  of  a 
particular  region  contain  the  same  mass. 
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TABLE  I:  Total  power  output  and  pump  rates  of  important  photoresonant 

transitions  in  Ne-like  krypton  I(j,fjo)T  notation].  The  mass  of  the  source 

I  ^  J 

plasma  is  273  ug/cm,  the  temperature  is  500  eV,  and  the  ion  density  is  5  x  10 

_3 

cm  .  The  4d  level  is  representative  of  a  number  of  fine  structure  split  levels 
that  are  lumped  according  to  statistical  equilibrium. 


Transition 
to  (ls^2p^) 
from 

Energy 

(eV) 

Oscillator 

Strength 

Power 

(lO^^W/cm-eV) 

Pump  Rate 
P(10^  sec~^) 

3s  (3/2, 1/2)^ 

1652 

0.136 

1.72 

9.67 

3s  (1/2, 1/2)^ 

1708 

0.095 

1.72 

6.78 

3d  (3/2,5/2)j 

1790 

1.66 

1.18 

74.5 

3d  (1/2, 3/2)^ 

1846 

2.0 

1.33 

98.6 

4d 

2309 

0.88 

0.32 

8.29 
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2p54d 


Fig.  1  Energy  level  diagram  notation)  and  oscillator  strengths 

for  important  lasing  and  resonance  transitions  for  neonlike  krypton. 
The  4d  level  is  representative  of  a  number  of  fine  split  levels  that 
are  lumped  according  to  statistical  equilibrium. 
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keV-sec 


gs/(cm-keV-sec 


1022 
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lO'a 
li  lO'S 
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1.6400  1.6450  1.6500  1.6550  1.6600 

Energy  (keV) 

Fig.  4  Emitted  spectral  profile  of  3s(3/2,l/2)j  -  ground  state  transit! 
of  neonlike  krypton  emitted  from  the  pump  plasma. 
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Effect i ve  Charge 


Fig.  5  The  effective  charge  as  a  function  of  density  and  temperature  oi 
lasant  plasma.  Radiation  from  the  pump  plasma  is  turned  off. 


5: 


3s  (J=0-1 .  172  A)  Gain 


Fig.  7  3p-3s  (J=0-1)  gain  coefficients  as  a  fun^-tion  of  temperature  for  the 

inner  zone  of  a  10^^  cm  ^  ion  density  lasant  plasma.  Radiation  from 
the  pump  plasma  is  turned  off.  The  radius  of  the  lasant  plasma  is 


0.1  cm 
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3p-3s  (J=2-1.  186  A)  Gain  (cnr 


Temperature  (eV) 

0 

Fig.  8  3p-3s  (J=2-l)  186  A  gain  coefficients  as  a  function  of  temperature 

and  density  of  the  lasant  plasma.  The  pump  plasma  is  present  but 
there  is  no  attenuation  of  the  Ne-like  ground  state  pump  rates  to 
the  3s,  3d  and  4d  levels.  The  radius  of  the  lasant  plasma  is  0.1  cm. 
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Fig.  9  3p-3s  (J=2-l)  186  A  gain  coefficients  as  a  function  of  temperature 
and  density  for  the  outer  zone  of  the  lasant  plasma.  The  pump 
plasma  is  present  and  there  is  attenuation  of  the  Ne-like  ground 
state  pump  rates  to  the  3s,  3d,  and  4d  levels.  The  radius  of  the 


lasant  plasma  is  0.1  cm. 
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Fig.  10  3p-3s  (J=2-l)  186  A  gain  coefficients  as  a  function  of  temperature 
for  the  outer  zone  of  the  low  density  lasant  plasma.  The  pump 


plasma  is  present  and  there  is  attenuation  of  the  Ne-like  ground 
state  pump  rates  to  the  3s,  3d,  and  4d  levels.  The  radius  of  the 
lasant  plasma  is  0.1  cm. 
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Ultra  Short  Vulaa  Xiaaar  Producad  Pluorlna  Plaaaa 

J.  Davis,  R.W.  Clark,  and  J.  Giuliani 
Radiation  Hydrodynamics  Branch 
Plasma  Physics  Division 
Naval  Research  Laboratory 
Washington,  D.C.  20375 

Abstract 

The  interaction  of  an  intense  ultra  bright  laser  pulse  with  a  planar 
fluorine  target  is  investigated.  The  interaction  is  numerically  simulated 
using  a  1-D  time-dependent  non-LTE  radiation  hydrodynamics  model,  including 
the  effects  of  the  ponderomotive  force,  self-consistently  coupled  to  a  wave 
equation  describing  the  beam-plasma  coupling.  The  target  is  shown  to  be 
highly  reflective  for  the  intensities  considered  but  still  allows 
sufficient  energy  to  be  coupled  to  produce  substantial  heating.  A  time 
history  of  the  evolving  plasma  is  presented  and  discussed  as  well  as  the 
role  played  by  the  ponderomotive  force.  The  feasibility  of  creating 
conditions  conducive  to  population  inversions  and  gain  in  selected 
transitions  in  helium-  and  hydrogen-like  fluorine  is  also  investigated. 

Introduction 

The  interaction  of  nanosecond  pulsed  lasers  with  a  variety  of  targets 
has  been  under  intensive  investigation  for  over  three  decades.  The  major 
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emphasis  of  the  research  was  focused  on  laser  fusion  experiments  dedicated 
to  creating  a  laser  produced  plasma  with  self  sustained  fusion  reactions. 
Although  interest  remains  high  for  the  eventual  success  of  the  laser  fusion 
program,  a  number  of  other  activities  are  emerging  as  a  result  of  recent 
technological  advances  in  laser  development.  One  of  these  activities  is  due 
to  the  application  of  chirped  pulse  amplification  making  it  possible  to 
produce  ultra  short  pulse  lasers  with  pulsewidths  of  10' s  to  100' s  of 
femtoseconds  and  intensities  in  excess  of  10^®  watts/cm^.  This  newly 
evolving  technology  is  introducing  a  number  of  innovative,  interesting  and 
exciting  activities  involving  the  interaction  of  ultra  short  pulse  lasers 
with  matter,  including  x-ray  laser  research,  laser  produced  hohlraum 
physics,  harmonic  generation,  subnanosecond  spectroscopy,  and  the 
interaction  of  intense  laser  fields  with  individual  atoms  and  molecules, 
just  to  mention  a  few.  The  interaction  of  these  intense  lasers  with  a  planar 
target  makes  it  feasible  to  investigate  solid  state  plasmas  in  the  strongly 
coupled  plasma  regime  and  preliminary  results  suggest  the  possibility  of 
tailoring  the  x-ray  pulse  emitted  from  the  hot  dense  plasma  before 
hydroi^  'iiamic  blowoff  begins  to  disassemble  the  target.  Once  expansion  of 
the  high  density  blowoff  plasma  occurs  it  can  cause  a  rapid  cooling  of  the 
leading  edge  of  the  plasma  resulting  predominantly  in  collisional 
recombination  into  the  upper  excited  levels  of  the  recombined  ion  creating 
the  possibility  of  population  inversions.  In  this  paper  we  investigate  the 
interaction  of  an  ultra  short  laser  pulse  with  a  planar  fluorine  target 
with  the  focus  of  the  research  directed  toward  characterizing  the  radiative 
properties  of  the  plasma  as  well  as  determining  whether  the  plasma  can 
support  population  inversions  and  lasing  in  selected  transitions  in  highly 
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ionized  excited  states  as  a  result  of  recombination  result  "ig  from  rapid 


cooling. 


Model 

The  numerical  simulations  are  based  on  a  1-D  time~dependent  non-LTE 
radiation  hydrodynamics  model  self-consistently  coupled  to  a  Helmholtz  wave 
equation  describing  the  interaction  of  the  laser  pulse  which  is  assumed  to 
be  normally  incident  on  the  target.  The  atomic  model  contains  all  the 
ground  states  and  selected  excited  states  distributed  throughout  various 
stages  of  ionization.  In  particular,  the  hydrogenlike  manifold  contains 
excited  states  up  to  and  including  the  n»10  level,  the  heliumlike 
ionization  stage  up  to  and  including  the  n-5  singlet  and  triplet  levels, 
and  the  lithiumlike  stage  contains  doublet  level  structure  up  to  and 
including  n«5.  The  remaining  ionization  stages  contain  only  ground  states. 
The  ionization  dynamics  is  based  on  a  time -dependent  collisional  radiative 
model  self-consistently  coupled  to  a  probability  of  escape  radiation 
transport  scheme.  Ionization  lowering  is  included  to  account  for  the 
disappearance  of  excited  levels  in  the  very  dense  regions  of  the  plasma. 
The  spectral  line  profile  functions  are  represented  by  Voigt  functions.  For 
an  in-depth  discussion  of  the  non-LTE  radiation  hydrodynamics  model  the 
reader  is  referred  to  references  1  and  2.  The  absorption  of  a  laser  pulse 
normally  incident  on  the  target  is  calculated  by  solving  the  Helmholtz  wave 
equation  for  a  multi-layered  medium.  In  this  way  we  self-consistently 
characterize  the  incident,  reflected  and  transmitted  waves  in  each 
numerical  zone  in  the  medium.  The  conductivity  was  calculated  using  the 
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Drude  model.  The  initial  ionization  was  produced  by  multiphoton  effects  but 
was  rapidly  dominated  by  inverse  bremsstrahlung  heating. 


Results 


The  reflective  properties  of  the  solid  density  planar  fluorine  target 
irradiated  by  a  normally  incident  248  nm  KrF  laser  pulse  with  a  peak 
intensity  of  3X10^"^  watts/cm^  and  a  pulse  duration  of  650  femtoseconds  is 
shown  in  Fig.l  as  a  function  of  time.  The  incident  laser  intensity  I  is 
normalized  to  Iq  (  in  this  case  lQ“3xlO^'^watts/cm^ )  .  Most  of  the  incident 
laser  energy  is  reflected  from  the  surface;  about  97%  as  shown  by  the  curve 
labeled  Reflect.  The  normalized  time  integrated  energy  absorbed  is 
represented  as  the  (  J (Incident-Reflected) dt/ J (Incident) dt  )  and  reaches  a 
peak  value  of  about  3  %  at  1/2  psec.  In  Fig. 2  the  energy  deposited  into  the 
target  at  0.755  psec  (  which  is  at  the  end  of  the  laser  pulse  )  is  shown  as 
a  function  of  position  (  the  original  position  of  the  target's  surface  was 
0.5)/)  .  Also  shown  is  the  ratio  averaged  over  an  optical  cycle  of  the  local 
electric  field,  <E>,  to  the  amplitude  of  the  incident  field,  Eq.  The 
oscillatory  nature  of  <E>  simply  mirrors  the  alteration  of  the  wavenumber  k 
due  to  target  absorption  -  reflection.  The  behavior  of  the  mass  density,  p, 
is  also  shown  as  is  ^critic'  ''hich  where  the  plasma  frequency  equals  the 
laser  frequency.  Note  the  increase  in  density  at  about  0.3)/  due  to  the 
shock  wave  launched  into  the  material  as  a  result  of  the  de^..jition.  The 
laser  energy  is  not  completely  reflected  at  P-j-ritic  penetrates  a  skin 
depth  further  into  the  target . 

xne  competition  between  the  gas  and  ponderomotive  pressure  is  shown 
in  Fig. 3  near  the  peak  of  the  laser  pulse  at  0.355  psec.  The  ponderomotive 
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force  exhibits  behavior  directly  related  to  the  square  of  the  electric 
field  and  exceeds  the  gas  pressure  in  the  leading  edge  of  the  blowoff 
plasma.  The  sudden  increase  in  both  the  ponde r omot i ve  and  fluid 
accelerations  beyond  0.1  jt  is  due  to  the  low  density  in  the  region.  At  this 
intensity  the  ponderomotive  force  plays  a  substantial  role  in  the  plasma 
evolution.  Figs.  4  and  5  show  the  dynamic  evolution  of  the  plasma  at  two 
times:  after  the  peak  of  the  laser  pulse  at  0.359  psec  and  the  end  of  the 
laser  pulse  at  0.759  psec,  respectively.  Note  the  change  in  the  horizontal 
axis  in  Figure  5  at  late  time.  The  break  in  the  velocity  curve  represents 
material  moving  outward  (  righthand  portion  )  and  inward  (  lefthand  portion 
) ,  respectively  .  The  peak  temperature  occurs  around  the  peak  of  the  laser 
pulse  and  is  in  excess  of  2  keV.  At  later  times  the  plasma  is  cooler,  and 
shows  a  pronounced  drop  in  temperature  near  1.0^.  The  dE/dtj-^^j  curve 
represents  the  total  radiated  power  density  over  all  wavelengths,  but  is 
certainly  dominated  by  short  wavelengths  where  the  plasma  is  hot . 

Emission  spectra  are  shown  in  figures  6  and  7  for  two  times:  after 
the  peak  of  the  laser  pulse  at  .359  psec  and  at  the  end  of  the  pulse  at 
.759  psec,  respectively.  The  earlier  spectrum  (fig.  6)  is  hotter,  and 
there  is  less  K-shell  recombination  radiation.  Near  the  peak  of  the  laser 
pulse,  about  20%  of  the  radiated  power  comes  from  spectral  lines 
(predominately  from  the  resonance  lines  of  hydrogen-  and  helium-like 
fluorine) ,  about  40%  comes  from  recombination  (largely  K-shell)  and  about 
40%  from  brems st r ahl ung ,  Most  of  this  radiation  emanates  from  the 
"transition"  region,  where  the  temperature  is  several  hundreds  of  eV,  and 
the  plasma  is  near  solid  density.  Because  this  region  is  of  nearly 
infinitesimal  thickness,  radiation  moving  in  the  forward  direction  is 
nearly  optically  thin.  The  latter  spectrum  (fig.  7)  is  characteristic  of  a 
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cooler,  denser  plasma.  However,  it  is  still  dominated  by  hydrogen-  and 
helinm-like  radiation,  and  the  fraction  of  the  emission  from  spectral 
lines,  recombination  and  bremsstrahlung  is  approximately  the  same. 
Finally,  the  time-integrated  emission  spectrum  is  shown  in  figure  8.  Not 
surprisingly,  it  is  very  similar  to  the  spectrum  near  the  peak  laser 
deposition  (fig.  6) . 

Finally,  for  the  conditions  considered  here  population  inversions 
were  not  observed  in  any  of  excited  states  of  hydrogen-,  helium-,  or 
lithium-like  fluorine.  It  may  be  possible  to  create  population  inversions 
leading  to  gain  at  higher  laser  intensities  and/or  different  geometries 
with  higher  Z  materials.  Both  possibilities  will  be  investigated. 
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Fig.  1.  Time  history  of  the  laser  deposition. 

Incident  laser  intensity  Itaier  is  normalized  with  respect 
to  the  peak  intensity  /<,.  The  reflection  coefficient  Reflect 
is  the  reflected  intensity  Ir^/  divided  by  the  incident  Ihater- 
The  normalized  time  integrated  energy  absorbed  is  given  by 
Ji^Later  -  I  Re  J  J Laterdi  .  It  is  plotted  xlO  for  legibifity. 
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Fig.  2.  Laser  electric  field  at  0.755  picoseconds. 

The  local  electric  field,  averaged  over  an  optical  cycle,  <  E  > 
is  shown  normalized  to  the  amplitude  of  the  incident  field  Eo- 
The  mass  density  p  is  shown  in  a  stepwise  manner  to  indicate 
the  spatial  resolution  of  the  simulation  in  the  deposition  region. 
The  logarithm  of  the  rate  of  energy  deposition  into  the  target 
Edep  is  plotted,  and  the  location  of  the  critical  surface  Rcrit 
is  shown. 
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Fig.  3.  Ponderomotive  acceleration  at  0.355  picoseconds. 

The  magnitude  of  the  local  ponderomotive  force 
is  compared  with  the  magnitude  of  the  fluid  (pressure  gradient) 
forces  just  after  the  peak  of  the  laser  pulse.  The  corresponding 
fluid  accelerations  are  plotted. 
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Fig.  4.  Hydrodynamic  profiles  at  0.359  picoseconds. 

Spatial  profiles  of  mass  density,  temperature,  pressure,  velocity 
and  radiative  power  are  plotted  near  the  peak  of  the  laser  pulse 
These  quantities  are  shown  on  a  logarithmic  scale,  normalized 
as  indicated  above. 
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Fig.  5.  Hydrodynamic  profiles  at  0.759  picoseconds. 

Spatial  profiles  of  mass  density,  temperature,  pressure,  velocity 
and  radiative  power  are  plotted  near  the  end  of  the  laser  pulse. 
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Fig.  6.  Emission  spectrum  at  0.359  picoseconds. 

Instantaneous  emission  spectrum  after  the  peak  of  the  laser  pulse. 
Hydrogen-like  Ly-a  and  Ly-/?  and  hehum-like  He-o  transitions  are 
indicated. 
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Fig.  8.  Time-integrated  emission  spectrum. 

Emission  spectrum  integrated  through  0.959  psec.  Note  the  simil¬ 
arity  to  the  spectrum  near  the  peak  of  the  laser  pulse. 
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Model  for  Time  Dependent  Radiative  Transfer  in 
Short  -  Pulse  Laser  Heated  Plasmas 

J.  P.  Apruzese 

Radiation  Hydrodynamics  Branch 


Abstract 


Multiple  scattering  of  line  photons  can  significantly  broaden  the  x-ray 
pulse  emitted  from  ultrashort  pulse  laser-heated  plasmas.  Time-dependent 
radiation  transport  is  a  potentially  important  effect  in  determining  both 
pulsevidth  and  intensity.  A  numerical,  fully  time  dependent  radiation 
transport  model  has  been  developed  and  benchmarked  using  the  two  level  atom 
formalism.  It  is  shown  that  two  Gaussian  quadrature  rays  provide  sufficient 
accuracy,  greatly  enhancing  its  computational  efficiency. 
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INTRODUCTION 


The  advent  and  rapid  development  of  powerful  subpicosecond  lasers  has 

stimulated  a  new  and  very  recent  generation  of  experiments  producing 

extremely  short  lived  plasmas^”^^  (for  a  review  see,  e.g.,  Ref.l).  Numerous 

uses  of  such  experiments  are  envisioned,  including  creation  of  ultrashort  x- 

ray  pulses, pumping  of  x-ray  lasers  by  field  ionization, x-ray 

17  18 

inner  shell  photoicnization,  ’  rapid  recombination  or  resonant 

photoionization.  Subpicosecond  x-ray  pulses  would  enable  unique  time  - 

resolvtfd  x-ray  imaging  and  scattering  studies  of  plasmas,  solids,  and 

biological  subjects.  Additionally,  production  of  coherent  XUV  radiation  by 

19 

backcattered  harmonics  from  powerful  short-pulse  lasers  has  been  proposed. 
Such  ultraintense  laser  radiation  will  also  facilitate  basic  investigations 
and  c'.scovery  of  fundamental  atomic  processes  in  the  presence  of  field 
strengths  equaling  or  exceeding  those  due  to  the  Coulomb  field  of  the 
nucleus. 

As  expected  in  such  a  burgeoning  field  in  its  relative  infancy, 
contradictions,  questions,  and  apparent  tradeoffs  have  emerged  in  the  initial 
series  of  experiments  carried  out  at  several  laboratories.  For  instance,  it 
is  still  far  from  settled  whether  the  ma’n  i^curce  of  x-rays  in  such  highly 
nonequilibrium  plasmas  will  be  lines  or  co.itinuun,  and  the  atomic  number  of 
the  irradiated  plasma  can  be  expected  to  play  a  major  role  in  determining  the 
principal  radiative  processes.  One  interesting  phenomenon  which  has  been 
both  identified  and  quantified  is  the  effect  of  a  laser  prepulse  on  the 
characteristics  of  the  emitted  x-rays.  Cobble  et  al.,^  and,  more  recently, 

9 

Teubner  and  co-workers  report  that  the  presence  of  a  laser  prepulse  has  a 
strongly  positive  effect  on  the  plasma  x-ray  yield  and  hence  conversion 
efficiency.  This  occurs  whether  the  prepulse  is  of  nsec,^  or  sub-psec^ 
duration.  Of  course,  inherent  in  the  prepulse  is  that  the  preformed  plasma 
has  time  to  expand.  As  pointed  out  in  Refs.  13-15,  hydrodynamic  expansion 
and  thermal  conduction  are  two  of  the  most  potentially  important  cooling 
mechanisms  in  such  plasmas.  The  greater  the  plasma  scale  length,  the  more 
these  mechanisms  tend  to  broaden  the  x-ray  pulse.  A  tradeoff  between  yield 
and  pulsewidth  thus  appears  to  be  indicated  by  recent  data.  Another  process 
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which  can  effect  both  the  intensity  and  the  width  of  the  x-ray  pulse  is 

radiative  transfer.  As  pointed  out  in  Ref.  15,  as  line  radiation  becomes 

more  important,  the  potential  for  significant  opacity  effects  increases. 

This  is  not  expected  to  be  important  for  submission-sized  plasmas  dominated 

13 

by  continuum  radiation.  However,  at  solid  density,  strong  resonance  lines 
can  be  very  optically  thick  even  for  such  small  path  lengths.  At  scale 
lengths  of  microns  or  longer,  which  have  been  experimentally  demonstrated  to 
produce  the  best  laser-to  x-ray  conversion  efficiency,  opacity  is  even  more 
likely  to  play  a  significant  role. 

In  the  present  work,  the  effects  of  opacity  on  the  principal  resonance 

line  of  He-like  A1  (a  commonly  employed  target)^’^’^  are  examined  in  detail 

using  a  fully  time-uependent  radiative  transfer  model  accounting  for  the 

effects  of  photon  transit  time  on  the  pulsewidth.  A  two-level-atom  formalism 

is  initially  adopted  to  benchmark  the  results.  The  properties  of  the  single 
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line  have  been  chosen  to  mimic  the  characteristics  '  thie  He-like  A1  Is  - 

ls2p  Pj,  line  at  7.757  A.  As  mentioned  above,  A1  (Z=13)  has  been  a 

frequently  used  target  in  ultrashort  pulse  laser  heating  experiments.  The 

elements  of  adjacent  atomic  number.  Mg  (2=12)  and  Si  (Z=1A)  have  also  been 
3  7  8 

popular  choices.  ’  ’  .  Therefore,  the  choice  to  model  in  detail  the 
strongest  A1  K  shell  line  is  clearly  of  experimental  relevance. 

RADIATION  TRANSPORT  MODEL 


The  plasma  conditions  considered  in  the  present  work  are  such  that  only 
spectral  lines  are  likely  to  be  optically  thick.  The  transport  of  radiation 
is  accounted  for  by  the  fully  time  dependent  equation  of  transfer  at 
frequency  v,  along  a  ray  whose  coordinate  is  denoted  by  x. 
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In  Eq.  (1)  I^  is  the  specific  intensity,  and  k^  are  the  emission  and 
absorption  coefficients,  respectively,  and  c  is  the  velocity  of  light.  For  a 
derivation  of  this  equation  see,  e.g..  Ref.  22.  By  convention  in  Eq.(l)  the 
radiation  is  traveling  in  the  direction  of  increasing  x.  Consider  the  space, 
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time  point  (x^,t^).  Looking  backwards  along  the  ray  x,  the  quantities 

determining  I^(x^,t^)  are  and  k^  encountered  at  retarded  times  t^, 

x)  =  t^  -  (x  -x)/c.  The  optical  depth  between  x  and  x  at  time  t  is 
o  o  o  o 


■  o 

t(x  ,t  ,x)  =1  k  (t  )  dx 
o  o’  J  V  '  r 

x 


(2) 


and  the  solution  to  Eq.  (1)  is 


’'b 

In  Eq.  (3),  Xj^  is  the  coordinate  of  the  plasma  boundary.  In  all  cases 
considered  in  the  present  work,  no  external  radiation  is  assumed  incident  on 
the  plasma. 

Given  the  emission  and  absorption  coefficients  (from  previous  time 
steps)  Eq.  (3)  is  solved  on  a  frequency  grid  of  31  points  per  transported  line 
and  a  spatial  grid  of  26-31  points.  In  the  numerical  solution,  and  k^  are 
linearly  interpolated  oetween  grid  points  in  space  and  time.  The  maximum  time 
step  employed  is  10  fsec.  To  compute  photoexcitation  and  escape  rates,  angle- 
averaged  intensity  and  net  flux  are  required,  necessitating,  in  general, 

the  use  of  more  than  one  ray.  In  the  present  work,  confined  to  plane-parallel 
geometry,  Gaussian  quadrature  is  employed  as  described  by  Chandresekhar  in 
Ref.  23.  Both  two  rays  (n=2,  m=A)  and  four  rays  (n=4,  m=8)  have  been  used 
comparatively  to  obtain  I^  and  as  described  in  the  (n,  m)  notation  in  Ref. 
23.  The  escape  rate  is  just  the  frequency  integral  of  F^.  The 
photoexcitation  rate  for  a  spectral  line  of  absorption  oscillator  strength  f 
from  lower  level  1  to  upper  level  u  is  given  by 
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In  Eq.(4)  uhe  integral  is  taken  over  the  line  profile  where  the 
normalization  is  given  by 
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Comparisons  of  solutions  given  by  the  2-ray  and  4-ray  Gaussian  quadrature 
approximations  appear  in  the  next  section. 


TWO-LEVEL  ATOM  SOLUTIONS 


To  aid  in  analysis  of  the  phenomena  revealed  by  fully  time-dependent 
coupled  transfer  and  kinetics  calculations,  and  also  as  a  means  of 
benchmarking  the  model,  it  is  most  appropriate  to  consider  first  the  case  of 
the  two-level-atom.  This  has  been  intensively  studied  for  the  time- 
independent  approximation  in  the  astrophysical  literature  (see,  e.g.,  Ref.  24 
and  references  therein).  For  reasons  described  in  the  Introduction  the 
properties  of  the  two-level-atom  considered  in  the  present  work  have  been 

chosen  to  correspond  to  those  of  the  principal  resonance  line  of  He-like  Al, 

2  1  ® 

Is  -ls2p  P. ,  whose  wavelength  is  7.757  A  and  spontaneous  decay  rate  is 

13  ^  -1 

2.78x10  sec.  .  Two  other  crucial  properties  of  the  line  are  a,  the 

broadening  parameter  for  its  assumed  Voigt  profile  and  e,  the  probability  that 

the  line  is  collisionally  quenched  by  electron  de-excitation  during 

21  -3 

reabsorption.  At  an  electron  density  of  10  cm  and  temperature  200  eV, 

these  parameters  are  obtained  as  a=0.05  and  c=0.029  from  multilevel  atomic 

data  for  Al  described  in  Ref.  25.  For  a  100  um  thick  plasma  slab,  the  optical 

depth  at  line  center  for  this  line  is  94,  measured  normally  across  the  full 

width.  These  assumptions  and  conditions  form  the  basic  framework  for  the 

initial  benchmark  calculations  to  be  described. 

Following  the  initial  collisional  excitation  of  the  upper  level  of  an 

optically  thick  line,  the  emitted  photon  can  be  absorbed,  usually  about  one 

optical  depth  away,  and  then  re-emitted  in  the  line  with  a  probability  (1-e) 

per  absorption.  The  probability  that  the  line  photon  escapes  the  plasma  (P^) 

rather  then  undergoing  absorption  can  be  calculated  by  averaging  over  the  line 

26  27  28 

profile  and  has  been  given  in  useful  analytic  forms  for  Doppler,  ’  Voigt, 

29  30 

and  Stark  ’  line  profiles.  The  probability  that  the  line  photon  remains  in 
the  plasma  following  emission  is  thus  (l-P^Xl-G).  The  photon  is  eventually 


eliminated  either  through  collisional  quenching  or  escape.  The  process  of 
absorption  and  re-emission  of  line  photons  is  generally  referred  to  as 
"scattering".  Evidently,  some  line  photons  will  escape  after  a  few 
scatterings  whereas  other  may  require  many  scatterings  or  end  up  being 
collisionally  quenched.  This  process  has  the  potential  to  significantly 
broaden  the  x-ray  pulse  emitted  from  short-pulse  laser-heated  plasmas. 

To  examine  this  mechanism  in  detail  as  well  as  validate  the  numerical 
model  the  following  initial  excitation  conditions  (t=0)  are  assumed.  The  100 
ym  thick  plasma  modeling  only  the  A1  XII  principal  resonance  line  is  initially 
entirely  in  the  ground  state  except  for  the  innermost  12X  which  contains  a  1% 
population  of  upper  level  excited  states,  as  illustrated  in  Fig.  1.  Very 
approximately,  this  corresponds  to  the  laser  heating  of  a  preformed  plasma 
near  the  critical  surface.  More  importantly,  it  provides  a  well-defined 
escape  probability  as  the  photons  have  a  common  spatial  origin.  This  permits 
analytical  comparison  with  the  numerical  results  and  facilitates  analysis  of 
the  physics.  Of  course,  even  when  one  employs  the  simplifying  two-level-atom 
approximation  there  is  more  to  the  calculation  than  evaluations  of  the 
photoexcitation  rates  as  per  Eqs.  (1-5).  The  full  equation  for  the  upper 
level  population  is 


dN 
_ u 

dt 


"e  -  "u  ‘V*  %1> 


(6) 


In  Eq.  (6)  W  is  the  photoexcitation  rate  as  given  in  Eq.  (4),  A  is  the 
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spontaneous  decay  rate  (2.78x10  sec  )  and  the  C's  refer  to  electron 
collisional  excitation  and  de-excitation  of  the  upper  level.  The  de¬ 
excitation  rate  is  easily  computed  from  the  condition 


e  =  0.029 


'ul 


(7) 
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ul 


The  excitation  rate  is  set  to  zero  to  simulate  the  effect  of  femtosecond- 

time-scale  plasma  cooling  given  the  intial  conditions  described  above,  and  fhe 

radiative  transfer  algorithm  for  calculating  the  W's,  Eq.  (6),  is  integrated 

31 

employing  the  method  of  Young  and  Boris. 

The  upper  level  fraction  as  a  function  of  space  and  time  is  illustrated 
in  Figs.  1  and  2  for  5  time  points  varying  from  0  to  7  psec.  These  numerical 
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solutions  to  Eqs  (1)  -  (6)  are  presented  for  both  4-ray  and  2-ray  Gaussian 

quadrature.  Note  first  the  excellent  agreement  between  the  n=2  and  n=4  cases, 

strongly  suggesting  that  2  rays  are  adequate.  Diffusion  of  the  upper  level 

population  through  the  plasma  through  line  photon  absorption  and  re-emission 

is  clearly  evident.  Of  considerable  interest  is  the  rate  of  decay  of  the 

upper  level  population  in  the  central  zone.  If  the  plasma  were  optically  thin 

this  population  would  decrease  with  an  e-fold  time  of  (2.78x10  )  =36  fsec. 

The  actual  e-fold  time  is  much  less  in  the  optically  thick  medium  because  of 

re-excitation  of  the  level  by  line  photons,  as  is  evident  in  Figs.  1  and  2. 

In  fact,  for  the  first  0.5  psec,  the  lifetime  of  the  upper  level  fraction  is 

0.7  psec,  and  this  lifetime  increases  to  1.1  psec  for  the  final  3  psec  covered 

by  Figs.  1  and  2.  A  frequently  used  analytic  approximation  for  line  transfer 

is  to  multiply  the  spontaneous  decay  rate  by  the  escape  probability, 

implicitly  assuming  that  "true"  de-excitation  occurs  only  a  fraction  of  the 

time.  According  to  the  calculations  of  Ref.  28,  the  escape  probability  for 

this  line  from  the  midpoint  of  the  plasma  is  8.9x10  .  Of  course,  collisional 

de-excitation  also  depletes  the  upper  level.  Employing  the  escape  probability 

12  -1 

dilation  method,  the  total  decay  rate  is  ~  A^j^  (P  +  e)  =  1.05x10  sec  , 
implying  a  lifetime  of  0.95  psec,  very  close  to  the  average  e-fold  time 
observed  in  Figs.  1  and  2.  Time  dependent  delays  in  re-excitation  result  in 
an  initially  more  rapid  population  falloff,  but  this  transient  phenomenon 
damps  out  over  7  ps,  resulting  in  good  correspondence  between  the  analytically 
estimated  and  numerically  computed  upper  level  decay  times. 

The  time  dependent  escape  of  the  line  photons  is  more  directly  relevant 
to  the  experimental  achievement  of  optimally  short  plasma  x-ray  pulses.  For 
the  two  level  atom,  this  quantity  may  also  be  analytically  estimated  and 
compared  with  the  present  numerical  calculations.  Ignoring  the  light  travel 
time  between  scatterings,  each  absorption  and  re-emission  requires  a  time 
equal  to  the  lifetime  L  of  the  upper  level  (36  fsec  in  this  case).  The 
probability  that  a  photon  escapes  after  an  excitation  is  (1-c)  P  ;  the 
probability  of  escape  by  the  N*’  scattering,  is  the  sum  of  the 
probabilities  of  escape  after  initial  excitation  and  each  of  the  subsequent 
scatterings,  viz. 


N, 

c  N+n 

P  =  (l-e)P 
n  '  e 

1+E  (l-E)"(l-P^)" 

.  n=l  ®  . 

=  (l-E)Pg 

l-S 

l-S 

(8) 


83 


where  in  Eq.(8)  S=(l-e)(l-Pg) ,  the  probability  that  the  photon  remains  in  the 
plasma  after  an  excitation.  The  time  required  for  N  scatterings,  i.e., 
associated  with  fractional  escape  P^j,  is  simply  the  level  lifetime,  36  fsec, 
multiplied  by  N.  In  the  case  under  consideration,  e=0.029,  P^=8.9  x  10,”^  and 
the  analytic  results  obtained  from  Eq.  (8)  are  plotted  along  with  the 
numerical  solutions  in  Fig.  3.  As  indicated  in  that  figure,  one  of  the 
plotted  numerical  solutions  ignored  the  time  dependent  term  in  Eq.  (1), 
equivalent  to  assuming  infinite  light  velocity.  As  expected,  it  is  this 
numerical  curve  which  most  closely  approximates  the  analytic  solution  (which 
also  ignores  the  light  travel  time  between  scatterings).  Note  that  the 
numerical  solutions  including  full  time  dependence  were  obtained  for  2  and  4 
Gaussian  rays,  and  the  2-ray  solutions  are  again  in  excellent  agreement  with 
those  of  4-rays  which  are  of  course  twice  as  computationally  expensive.  It  is 
clear  that  the  economical  2-ray  quadrature  is  adequate. 

The  asymptotic  or  ultimate  escape  probability  for  a  line  photon  is 
obtained  by  taking  the  limit  of  Eq.  (8)  as  N  ->  ®.  For  the  presently 
considered  conditions  a  value  of  0.23  is  obtained,  as  also  seen  in  the 
analytic  curve  of  Fig.  3.  Most  photons  (77Z)  do  not  escape,  since  the 
quenching  probability  exceeds  that  of  escape  by  a  factor  of  -  3.  The 
numerical  solutions  exhibit  two  important  departures  from  the  analytic  result. 
First  the  photons  do  not  emerge  immediately  since  170  fsec  is  required  to 
traverse  the  50  wm  to  the  outer  edge  of  the  plasma.  Second,  about  28%  of  the 
initially  created  line  photons  escape  rather  than  the  23%  estimated 
analytically.  This  is  not  particularly  surprising  since  the  photons  diffuse 
toward  the  boundary  spatially  (as  seen  in  Figs.  1-2)  and  therefore  do  not 
strictly  obey  an  escape  probability  calculated  from  the  midplane.  It  is  also 
seen  from  Fig. 3  that  the  pulsevidth  obtained  from  the  full  numerical  solutions 
which  include  the  time-dependent  term  exceeds  that  yielded  by  dropping  the  1/c 
(3I/3t)  term.  For  the  time-dependent  case,  3.1  psec  is  required  for  90%  of  the 
photons  to  escape  whereas  2.4  psec  is  needed  when  this  effect  is  ignored. 
Clearly,  the  light  travel  time  can  affect  the  x-ray  emissions  when  very  short 
pulsewidths  are  considered  and  multiple  photon  scattering  is  significant  in 
the  overall  x-ray  energy  budget.  This  effect  can  be  examined  in  more  detail 
by  coupling  the  time-dependent  radiative  transfer  to  a  multilevel  ion  kinetics 
model . 
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FIG.  2  As  in  Fig.  1,  except  that  times  of  A.O  and  7.0  psec  after  initial 
excitation  are  illustrated. 
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Abstract 


We  present  CRE  gain  calculations  for  the  important  J=2-l  and  J=0-1  3s-3p 
transitions  of  neon-like  selenium  using  state-of-the-art  excited  state 
couplings  and  self-consistently  calculated  ionization  and  excited  state 
abundances.  An  atomic  model  was  put  together  containing  all  the  most 
significant  processes  impacting  the  gain  calculation  of  these  transitions  in 
order  to  determine  the  sensitivity  of  the  gain  ratio  of  the  J=0-1  to  the  J=2-l 
lines  to  the  monopole  excitation  cross  sections  in  CRE.  It  is  seen  that 
variations  of  a  factor  of  2  or  3  of  these  cross  sections  are  all  that  are 
needed  to  reverse  the  calculated  J=0-l/J=0-2  gain  ratios.  For  example,  a 
decrease  in  the  direct  excitation  rate  from  the  ground  to  the  (1/2, 1/2)^,  3p 
level  by  a  factor  of  three  is  sufficient  to  reverse  the  gain  ratio  for  the 
J=0-1  line  at  182  A°  to  the  J=2-l  line  at  210  A°  and  move  it  into  agreement 
with  the  experimentally  observed  ratio  for  these  transitions.  Experimental 
observations  would  suggest  that  for  low  Z  elements  (Cu,Zn),  this  ratio  is 
larger  than  one,  while  for  high  Z  elements  (Mo, Sr),  it  is  less  than  one.  At 
issue  is  where  the  cross-over  occurs.  We  draw  attention  to  another  instance 
where  the  discrepancy  between  theoretically  calculated  behavior  and 
experimental  behavior  can  be  attributed  to  the  large  size  of  the  calculated 
monopole  excitation  rate.  We  also  present  the  calculated  temperature  and 
density  behavior  of  the  J=0-1  and  J=2-l  gains  and  discuss  how  this  behavior 
might  be  used  to  determine  the  adequacy  of  present  atomic  models  to  explain 
experimental  neon-like  X-ray  laser  observations. 
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I.  INTRODUCTION 


1  2 

Since  the  first  demonstrations  of  soft  X-ray  lasing  in  neon-like  selenium  ’  , 

a  number  of  experimental  and  theoretical  developments  have  occured  leading  to 

improved  gain  calculations  and  measurements  and  to  the  achievement  of  even 

shorter  wavelength  neon-like  x-ray  lasers.  The  population  inversion  scheme 

employing  neon-like  selenium  has  been  thoroughly  investigated  because  the 

energy  multiplet  structure  is  relatively  simple,  and  neon-like  selenium  X-ray 

lasers  have  been  carefully  studied  experimentally.  However,  several  important 

issues  regarding  the  interpretation  of  experimental  results  by  these 

calculations  remain  to  be  adequately  resolved.  The  3s-3p  transitions  in  the 

neon-like  system  were  considered  a  good  source  for  lasing  because  the  2s  2p  3s 

lower  lasing  levels  can  decay  to  the  ground  state  by  a  dipole  allowed 

radiative  transition  in  less  than  a  picosecond  whereas  the  upper  metastable 
2  5 

2s  2p  3p  lasing  levels  can  do  so  much  more  slowly  only  by  a  quadrupole 

transition  and  by  collisional  deexcitation.  Neon-like  selenium  was  considered 

to  be  a  very  good  initial  candidate  ion  for  experimental  study  since  the  first 

theoretically  predicted  gains  for  the  plasma  conditions  achievable  in  the 

2 

early  experiments  seemed  to  peak  around  selenium  .  In  recent  years  theoretical 

calculations  and  experimental  measurements  have  been  carried  out  for  many 

different  Ne-like  ions.  In  order  to  calculate  the  population  inversion  between 
2  5  2  5 

the  2s  2p  3p  and  2s  2p  3s  excited  states  all  possible  mechanisms  responsible 
for  populating  and  depopulating  these  states  have  been  considered.  The 
different  theoretical  models  that  have  been  developed  have  varying  strengths 
and  weaknesses  because  they  include  detailed  calculational  data  for  some  of 
the  atomic  processes  while  other  processes  are  incorporated  only  in  an 
approximate  way.  The  major  problem  addressed  by  recent  calculations  is  to 
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explain  the  experimental  findings  for  the  relative  strengths  of  the  strongest 
lasing  lines.  The  greatest  discrepancy  between  the  experimental  observations 
and  theoretical  calculations  that  still  puzzles  scientists  today  is  the  fact 
that  the  J^O-l  182  line,  which  is  calculated  to  have  the  highest  gain, 
showed  much  less  or  little  gain  when  measured.  Recent  improvement  in 

3 

experimental  procedures  and  diagnostics  and  more  complete  theoretical  models 

4 

and  calculations  utilizing  detailed  and  more  accurate  atomic  data  have 
lessened  this  discrepancy,  but  the  disagreement  persists.  In  accord  with  these 
developments,  the  work,  to  be  presented  in  this  paper  utilizes,  in  detail  and 
self-consistently,  an  atomic  model  that  includes  the  atomic  processes  that  are 
known  to  be  responsible  for  populating  or  depopulating  the  laser  levels  in 
order  to  analyze  and  understand  systematically  how  significantly  each  process 
affects  CRE  gain  calculations  of  the  neon-like  selenium  3s-3p  lines  in  order 
to  explain  why  the  J«0-1  discrepancy  may  not  be  such  a  mystery  at  all. 

There  are  two  neighboring  J-0  and  J=2  3p  levels  that  produce  important 
lasing  transitions  in  neon-like  Se.  Ve  will  focus  specifically  on  these  states 
to  quantify,  as  best  possible,  the  different  pumping  mechanisms  that  determine 
the  J»0-1  to  J«2-l  gain  ratios  and  the  sensitivity  of  these  ratios,  in 
particular,  to  the  monopole  excitation  rates  of  the  J=0  states  from  the  neon¬ 
like  ground  state.  Extensive  experimental  investigations  have  been  carried  out 
in  recent  years  with  the  goal  of  accurately  measuring  these  ratios  in  a 
variety  of  different  experiments.  These  investigations  include  not  only 
efforts  to  improve  upon  the  experimental  data  for  Se  but  other  Ne-like  ions  as 
well.  For  example,  experimental  data  on  Ne-like  yttrium^  was  obtained  from  the 
same  experiments^  in  which  Se  was  investigated.  However  interpretation  of  the 
measured  gains  was  ambiguous  because  of  line  coincidence  between  the  J=0-1  and 
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5 

one  of  the  J«2-l  transitions.  The  experimental  data  of  molybdenum  had 
similarity  with  the  Se  data  i.e,  the  J»2-l  lines  were  measured  to  have 
substantial  gain  whereas  the  J*0-1  line  analogous  to  the  182  A°  line  in  Se  did 
not  show  much  amplification.  In  their  observation  of  amplication  in  the  3s-3p 
lines  in  Ne-line  germanium,  researchers  at  the  Naval  Research  Laboratory 
(NRL)^  obtained  measurable  gain  for  the  J»0-1,  196  A°  line,  but  the  gain  for 
the  J*2-l  lines  were  higher  compared  to  this  J*0-1  line.  Thus  for  all  ions, 
the  relative  strengths  of  the  measured  gain  for  the  J=2-l  lines  were  larger 
and  sometimes  much  larger  than  the  measured  gain  of  the  J=0-1  lines.  However 
for  copper^,  these  ratios  were  reversed  (as  the  original  gain  calculations  for 
selenium  had  predicted).  The  J*0-1  line  at  221  A°,  analogous  to  the  182  A° 
line  in  Se,  showed  a  larger  gain  compared  to  the  two  J=2-l  lines.  More  recent 
experiments  on  zinc,  gallium,  arsenic,  and  selenium^  have  also  been  performed 
at  NRL.  Gain  measurements  could  not  be  obtained  for  Ne-like  gallium  because 
the  lines  were  too  weak.  There  were  also  no  observed  gain  in  arsenic  for  the 
J«0-1  line,  analogous  to  the  182  A°  line  in  Se.  In  the  same  experiment,  gains 
of  4.9  cm~^  were  measured  for  the  two  J=2-l  lines  for  Se  and  a  measurable  (but 
substantially  lower  than  theoretically  predicted)  gain  of  2.6  cm  ^  was 
reported  for  the  J«0-1,  182  A°  line.  In  the  same  experiment  however,  the  J=0-1 
line  for  zinc  showed  a  larger  gain  coefficient  than  the  gain  for  the  two  J=2-l 
lines.  These  observations  would  indicate  that  the  relative  strengths  of  the 
gains  for  the  J«2-l  are  indeed  smaller  compared  to  the  gain  of  the  J=0-1  line 
for  sufficiently  low  Z  ions. 

Recent  improved  experiments  in  Se  at  LLNL  achieved  a  measureable  gain  for 
the  J-0-1  182  A°  line.  The  results  of  this  experiment  are  in  much  closer 
agreement  with  the  new  calculations^  of  the  theoretical  group  at  LLNL.  The 
better  agreement  resulted  because  the  gain  of  the  J=2-l  lines  increased  and 
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the  gains  of  the  J>0-1  lines  decreased  in  their  calculation,  and  reanalysis  of 
the  gain  data  of  the  J^O-l  line  in  the  measurement  shoved  a  noteworthy  gain  of 
2.4  cm“^.  Very  recently  substantial  gain  has  been  measured  for  the  J=2-l 

O 

transition  of  Ne-like  strontium  .  The  J=0-1  transition,  analogous  to  the  182 
A°  line  in  Se,  was  not  observed  due  to  a  wavelength  overlap  with  a  Na-like  Sr 
line.  Thus,  from  all  tho  data  available  so  far  for  the  Ne-like  ions  there  is 
strong  indication  of  a  systematic  Z-scaling  behavior  of  the  gain  ratio  of 

Q 

J*0-1  to  J=2-l  lines.  This  scaling  was  first  noted  by  Whitten  et  al  .;  their 
calculated  gains  for  the  J=2-l  lines  increased;  whereas  the  J=0-1  gains 
decreased  as  a  function  of  Z.  The  difference  between  the  calculated  and 
measured  gain  ratios  appears  to  be  in  the  Z  value  at  which  the  ratio  becomes 
one.  In  order  to  understand  the  Z  scaling  of  the  gains  of  the  two 
different(J*2-l  and  J=0-1)  but  adjacent  transitions  and  particularly  the 
behavior  of  the  J-0-l/J«2-l  gain  ratio  as  a  function  of  Z,  one  has  to 
accurately  calculate  the  Z  dependence  of  the  relative  strengths  of  the 
different  pumping  mechanisms  that  predominantly  populate  the  upper  laser 
levels. 

Collisional  excitation  is  the  single  most  important  mechanism  for 
populating  the  J=0  levels  of  the  3p  multiplet.  The  evaluation  of  the  monopole 
excitation  rate  from  the  ground  2p  state  has  received  much  attention.  Thus 
most  atomic  models  include  a  well  calculated  collisional  excitation  rate  in 
their  gain  calculations.  Although  other  ionization  and  recombination  processes 
must  also  be  included,  their  influence  on  the  population  inversion  dynamics  is 
less  direct. 

It  was  also  realized  that  the  J=2  upper  lasing  states  could  be  pumped 
preferentially  by  dielectronic  recombination  (DR)  from  f lourine-like  ground 
state  relative  to  the  J^O  state  pumping.  Thus  DR  was  thought  to  have  a 


96 


significant  effect  on  the  gain  calculation  of  the  lasing  levels  especially  for 
non-equilibrium  plasmas,  and  the  inclusion  of  detailed  DR  rates  into  the 
kinetic  model  for  each  of  the  ni=3  multiplet  states  seemed  absolutely 
essential^^’^^.  However,  such  inclusion  still  did  not  fully  explain  the  "J=0" 
mystery  in  Se.  Early  calculations  contained  only  a  ground-to-ground  effective 
DR  rate.  The  theoretical  group  at  Livermore,  in  their  recent  improved 
calculation,  included  state  specific  DR  rates  to  each  of  the  lover  and  upper 
lasing  levels.  These  inclusions  and  other  improvements  in  their  atomic  model 
were  observed  to  be  responsible  for  the  better  agreement  between  their 

3 

calculated  and  observed  gains  . 

Another  dominant  mechanism  for  populating  the  3p  levels  is  deexcitation 
from  the  3d  excited  levels.  Like  DR,  3d  deexcitation  also  favors  the  pumping 
of  J«2  states  over  that  of  the  J«0  states.  It  was  pointed  out  by  Whitten  et. 
al.^^  that  collisional  and  radiative  cascade  from  the  excited  states  into  the 
3d  states  increase  the  gain  of  the  J=2-l  lines  by  15-20)(. 

12 

It  was  also  shown  that  the  rates  for  inner-shell  ionization  of  sodium¬ 
like  ions  could  be  comparable  to  electron  excitation  from  neon-like  ions  in  a 
nonequilibrium  plasma.  In  fact  inner-shell  ionization  also  pumps  3p  states 
while  favoring  the  J»2-l  inversion  over  the  J=0-1.  This  process  might  also 
figure  into  the  interpretation  of  experimental  findings,  especially  for  short 
pulse  experiments. 

In  this  work,  we  have  constructed  a  detailed  atomic  model  of  the  n=3 
states  of  neon-like  selenium  and  coupled  it  self-consistently  to  a  less 
detailed  but  fairly  complete  model  of  the  ioni'ation  dynamics  in  all 
ionization  stages  of  the  selenium  system.  The  model  makes  use  of  state-of-the- 
art  atomic  data  and  contains  all  of  the  processes  known  to  play  an  important 
role  in  the  excitation  dynamics  of  the  n=3  excited  states. 
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It  enables  us  to  carry  out  self-consistent  calculations  of  both  the 
ionization  abundances  of  the  neon-  and  fluorine-like  ground  states  and  of  the 
populations  of  each  of  the  26  multiplet  states.  Vith  it,  we  have  computed 
the  collisional  radiative  equilibrium  (CRE)  behavior  of  the  J>0-1  and  J32-1 
gain  coefficients  and  quantified  their  dependence  on  the  three  important 
competing  processes  for  populating  the  3p  states  from  the  ground  neon-  and 
fluorine-like  states  and  from  the  12  3d  multiplet  states.  Our  model  produces 
somewhat  different  relative  magnitudes  of  the  different  X-ray  laser 
transitions  in  neon-like  Se  than  other  models  that  have  been  reported  on,  and 
we  use  it  to  evaluate  the  relative  importance  of  the  different  pumping 
mechanisms  in  interpreting  the  gain  calculation  results. 

This  paper  is  structured  as  follows.  In  section  II,  the  atomic  model  that 
we  use  for  the  gain  calculations  is  described.  The  relative  strengths  of  the 
three  dominant  CRE  pumping  mechanisms  are  discussed  in  section  III.  These 
magnitudes  depend  not  only  on  the  relative  strength  of  the  collision  rates, 
but  on  the  relative  size  of  the  populations  being  pumped.  The  parametric 
dependence  of  the  J*0-1  and  J=2-l  gains  on  temperature  and  density  is  then 
presented  in  section  IV,  where  the  sensitivity  of  these  results  to  the 
strength  of  the  monopole  excitation  rates  is  discussed.  Finally,  in  section  V, 
we  summarize  the  results  of  these  calculations  and  discuss  how  they  bear  on 
the  problem  of  comparing  theoretical  gain  calculations  with  experimental 
results. 
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II.  ATOMIC  MODEL 


Our  atomic  model  was  developed  to  include  all  the  important  processes 

responsible  for  3s-3p  population  inversions  in  neon-like  selenium  but  in  as 

economical  yet  complete  a  fashion  as  possible.  Thus  to  compute  the  overall 

ionization  state  of  selenium,  we  began  with  the  hydrogenic  atomic  model 

described  in  Reference  13,  and  replaced  the  lumped  excited  state  structure  of 

the  hydrogenic  model  with  the  30  excited  states  shown  in  Fig.  1  in  order  to 

compute  the  detailed  excitation  dynamics  of  the  n=3  states  of  neon-like 
24+ 

selenium,  Se  .  Similarly,  the  fluorine-like  structure  of  the  hydrogenic 

model  was  replaced  by  a  lumped  state  version  of  Fig.  1.  As  we  can  see  from 

this  figure,  the  model  contains  all  the  multiplet  levels  (4  levels  of  3s,  10 

levels  of  3p  and  12  levels  of  3d)  of  valence  n=3  excited  states  of  immediate 

relevance  to  neon-like  lasing  as  well  as  the  3  inner  shell  excited  2s2p^3s, 

2s2p^3p,  2s2p^3d  states,  and  the  n>4  lumped  excited  state  of  indirect 

relevance.  Ve  include  all  the  n-3  multiplet  states  and  all  processes  coupling 

these  levels  explicitly  because  these  populations  support  large  inversion 

densities  only  when  the  levels  are  far  from  LIE,  i.e.,  only  when  the 

metastability  of  the  n^l  states  is  pronounced.  The  fluorine-like  ionization 

stage  was  modeled  using  tha  same  atomic  structure  as  the  neon-like  system  but 

in  a  less  detailed  fashion  since  only  its  ground  state  populations  play  an 

important  role  in  determining  the  n=3  neon-like  populations.  For  economy,  it 

contained  a  lumped  description  of  the  57  n=3  multiplet  states;  thus,  in  all, 

it  contained  only  eight  excited  states  :  the  An=0  excited  state,  the  3s,  3p, 

5  5  5 

and  3d  lumped  states,  and  the  2s2p  3s,  2s2p  3p,  2s2p  3d,  and  n:=4  states  of 
Fig.  1.  The  CRE  gain  calculations  consist  of  simultaneous  calculations  of  the 
populations  of  all  ionization  stages  surrounding  the  neon-like  ionization 
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stage  in  addition  to  the  populations  of  all  the  neon-like  excited  states  in 
Fig.  1. 

Figure  2  shows  the  lA  specific  J  levels  that  constitute  the  upper  3p  and 
lower  3s  multiplet  states  in  neon-like  selenium.  The  four  J=0-1  and  J=2-l 
radiative  decays  along  with  the  wavelengths  for  these  transitions  are  also 
shown  in  this  figure.  It  is  drawn  to  scale  and  defines  labels  for  each  of  the 
four  J*0-1  and  J*2-l  transitions,  which  are  the  focus  of  the  theoretical 
investigations  of  this  paper. 

The  processes  that  couples  these  levels  are  collisional  and  resonance 
excitation,  collisional  deexcitation,  radiative  decay,  radiative 
recombination,  and  dielectronic  recombination  both  from  the  fluorine-like 
ground  state  and  its  dn^O  excited  state.  It  is  crucial  to  obtain  these 
couplings  as  accurately  as  possible  in  order  to  calculate  the  population 
dynamics  of  these  states  accurately.  Ve  have  obtained  this  data  from  a  variety 
of  sources;  all  of  it,  to  the  best  of  our  knowledge,  is  state-of-the-art.  For 
example,  we  used  the  collision  strengths  calculated  by  Zhang  et  al.^^  to 
couple  the  n*3  states  to  the  neon-like  ground  state.  The  collisional 
excitation  rate  coefficients  were  then  obtained  by  fitting  the  collision 
strength  data  points  (see  Fig,  8)  and  using  the  fits  to  integrate  the  cross 
sections  over  a  Maxwellian  electron  distribution  for  each  required 
temperature.  The  resulting  rate  coefficients  to  the  most  important  upper  and 
lower  lasing  levels  are  shown  in  Fig.  3.  It  shows  by  how  much  the  collisional 
excitation  rates  to  the  upper  J=0  levels  are  larger  than  the  rates  to  the  J=2 
levels  or  to  the  lower  J=1  levels.  Collisional  ionization  and  recombination 
rates  from  the  3s  sublevels  were  obtained  from  Golden  et.  al.,^^  and  the 
ionization  and  recombination  rates  from  the  3p  and  3d  sublevels  were  taken 
from  Moore  et.  al.^^. 


100 


Because  collisional  deexcitation  from  the  3d  excited  levels  to  the  3p 

states  is  a  major  contributor  in  populating  the  J=2  upper  lasing  levels  and 

because  3p-3p  collisions  tend  to  equalize  the  gain  of  the  two  J=0-1  lines, 

it  is  important  that  these  collisional  couplings  be  accurately  calculated  in 

our  model.  We  have  used  the  relativistic  distorted  wave  3s-3s,  3s-3p,  3p-3p, 

3s-3d,  3p-3d,  and  3d-3d  collision  strengths  calculated  by  P.  Hagelstein  and  R. 
18 

Jung  .  Again,  ve  used  algebraic  fits  to  the  data  given  in  their  paper  to 

carry  out  the  integration  over  a  Maxwellian  electron  distribution  to  obtain 

the  rate  coefficients  coupling  the  n=3  states. 

Additional  coupling  of  the  neon-like  excited  states  to  the  fluorine-like 

system  occurs  by  radiative  and  dielectronic  recombination  (DR).  Only  a  total 

ground-to-ground  radiative  recombination  coupling  of  the  fluorine-like  to  the 

neon-like  ionization  stage  was  included  in  our  model.  The  DR  rates,  on  the 

19 

other  hand,  were  calculated  from  a  detailed  calculation  of  the  needed  doubly 

excited  energy  levels,  stabilizing  radiative  rates,  and  autoionization  rates 

using  the  HFR  (Hartree-Fock  with  relativistic  correction)  method  of  R.  D. 
20 

Cowan  .  The  inclusion  of  ground-to-ground  effective  DR  rates  into  the  rate 

equation  dynamics  had  been  argued  early  on  in  the  history  of  the  the  neon-like 

Se  x-ray  laser,  and  more  recently  the  importance  of  state  specific  (for  each 

excited  multiplet  level)  DR  processes  had  been  claimed  by  London  et.  al.  .  We 

have  taken  into  account  a  complete  set  of  doubly  excited  states  in  calculating 
19 

the  DR  rates  .  For  low  j.ying  doubly  excited  states,  the  DR  branching  ratios 
were  explicitly  calculated,  while  for  higher  Rydberg  states  the  states  were 

3 

lumped  and  the  1/n  falloff  of  the  DR  rates  was  used,  where  n  denotes  the 

principal  quantum  number  of  the  Rydberg  electron.  Our  total  ground  to  ground 

DR  rate  coefficients  compared  very  well  with  the  detailed  multi-configuration 

21 

Dirac-Fock  (MCDF)  calculation  M.  C.  Chen  .  The  DR  rate  coefficients  that  were 
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obtained  from  our  calculations  connecting  the  ground  state  of  F-line  selenium 
to  the  most  important  upper  3p  and  lover  3s  lasing  levels  are  shown  in  Fig.  2 
of  Ref.  19.  This  figure  shows  that  coupling  strengths  from  the  F-like  ground 
state  are  in  reverse  order  to  those  from  the  neon-like  ground  state,  i.e.,  the 
J=2  levels  of  the  3p  states  are  pumped  more  strongly  by  dielectronic 
recombination  (DR)  than  the  J*0  levels  are.  The  rate  coefficients  for  the 
missing  multiplet  level  2p^3p  (3/2, 3/3)^  in  that  figure  are  approximately  of 
the  same  magnitude  as  those  for  the  2p^3p  (1/2, 1/2)^  level. 

Besides  direct  collisional  excitation,  the  excitation  of  neon-like  ions 

can  proceed  via  intermediate  autoionizing  resonances  which  can  often  have 

contributions  comparable  to  or  even  larger  than  those  due  to  direct 

excitation.  These  resonances  are  produced  when  a  neon-like  ion  captures  a 

target  electron  into  a  doubly  excited  sodium-like  state,  which  may  either 

radiatively  decay  or  autoionize.  The  radiative  channel  leads  to  dielectronic 

recombination  to  sodium-like  states  but  when  doubly  excited  states  autoionize 

to  excited  neon-like  states,  they  act  as  intermediaries  in  the  resonance 

excitation  (RX)  of  the  neon  ion.  We  calculated  these  rates  in  the  resonance 

22 

approximation  following  Covan  in  a  way  that  was  very  similar  to  the  DR  rate 

calculation.  Their  importance  to  the  3s-3p  gain  calculation  is  illustrated  in 

Fig.  4,  which  shows  not  only  that  the  RX  rates  to  the  3s  levels  are  much 

larger  than  those  to  the  3p  levels  but,  by  comparison  to  Fig.  3,  that  the  3s 

RX  rates  are  larger  than  the  direct  excitation  to  these  states.  The  inclusion 

of  these  processes  therefore  will  increase  the  J=l,  3s  populations  slightly 

and  decrease  the  calculated  gain  coefficients.  When  we  compared  these  level 

23 

specific  RX  rates  with  those  of  Chen  ,  we  found  that  our  calculated  RX  rates 
to  the  3p  multiplet  levels  are  lower  than  his,  but  our  rates  to  the  3s 
multiplet  states  compare  quite  well  with  his  rates. 
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Because  accurate  DR  rates  are  also  important  to  the  calculation  of  the 
ionization  balance,  which  indirectly  has  a  large  influence  on  the  3s-3p  neon¬ 
like  gain  calculations,  ve  have  also  included  in  our  model  a  detailed 
calculation  of  the  DR  rates  from  the  0-like  ground  and  /ki=0  excited  states  to 
all  levels  of  the  singly  excited  F-like  ion.  The  inclusion  of  these  processes 
has  a  significant  impact  on  reducing  the  oxygen-like  fractional  populations, 
and  hence,  in  turn,  on  shifting  the  CRE  ionization  balance  and  influencing  the 
strength  with  which  the  neon-like  states  are  pumped. 

Finally,  to  complete  the  collisional  couplings  in  our  Ne-and  F-like 

models,  we  have  also  calculated  excitation  and  deexcitation  collision  rates 

between  all  the  other  lumped  states  of  our  model  and  each  of  the  lover 

multiplet  levels  by  again  integrating  cross-sections  over  Maxwellian  electron 

distributions.  The  corresponding  cross-sections  were  calculated  in  a  non- 

relativistic  distorted-wave  approximation  using  atomic  wave  functions  and  data 

20 

generated  by  the  computer  code  of  R.  D.  Covan  .  Intermediate  coupling  and 
configuration  mixing  were  taken  into  account.  The  rates  from  the  n=4  lumped 
state  to  all  lower  n=3  multiplet  states  were  obtained  from  individual  J-J' 
rates  with  the  assumption  that  the  population  of  individual  upper  J  levels  is 
proportional  to 
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III.  PUMPING  MECHANISMS 


Collisional  excitation  from  the  Ne-like  ground  state,  dielectronic 
recombination  from  the  F-like  ground  and  An=0  excited  states,  and  collisional 
deexcitation  from  the  3d  excited  states  are  the  most  important  competing 
atomic  processes  that  are  responsible  for  directly  populating  the  upper  3p 
lasing  levels.  Fig.  5  shows  an  energy  level  diagram,  which  applies  to  either 
of  the  two  neighboring  J=0  and  J=2  upper  lasing  level  pairs,  that  summarizes 
the  situation  in  which  these  three  pumping  mechanisms  compete  to  populate  the 
two  states.  The  dark  arrows  designate  which  of  the  two  rates  to  the  J=0  and 
J«2  states  is  the  strongest  in  each  case.  Note  that  the  J=0  level  is  pumped 
preferentially  only  by  collisional  excitation  from  the  -round  Ne-like  state, 
while  the  J=2  level,  on  the  other  hand,  is  pumped  preferentially  by  both 
collisional  deexcitation  from  the  3d  excited  state  and  by  recombination  (DR) 
from  the  F-like  ground  state.  Note  also  that  the  populations  of  the  3d  states, 
which  most  directly  affect  the  population  of  the  J=2  level,  are  themselves 
directly  populated  by  collisional  excitation  from  the  Ne-like  ground  state, 
and  by  direct  and  cascade  recombinations  as  indicated  in  this  figure.  Two  of 
the  3d  states  are  depopulated  by  strong  radiative  decays  to  the  ground  state, 
while  the  other  10  have  no  strong  decays  and,  like  the  3p  states,  are 
metastable  (see  Fig.  10). 

Fig.  5  partially  explains  why  the  gain  ratio  of  neighboring  J=0-1  to  J=2-l 
lines  can  turn  out  to  be  larger  or  smaller  than  unity.  The  relative  strengths 
of  the  two  gains  will  also  depend  on  the  relative  strength  of  the  three 
pumping  rates  into  the  J=0  and  J=2  levels.  Since  these  pumping  rates  are  each 
the  product  of  the  total  rate  coefficient  for  a  given  process  times  the 
population  of  the  state  undergoing  the  transition,  the  total  pumping  strength 
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depends  on  the  fractional  population  that  is  calculated  for  each  state  as  a 

function  of  the  plasma  conditions.  Thus,  on  the  one  hand,  the  total  pumping 

strengths  to  the  two  J*0  upper  levels  may  be  small  if  the  strong  collisional 

excitation  rate  coefficient  multiplies  a  small  fractional  population  of  the 

Ne-like  ground  state.  (This  effect  was  argued  early  in  the  history  of  the 

selenium  X-ray  laser^^’^^.)  Similarly,  the  total  pumping  strengths  to  the  J=2 

upper  levels  depend  not  only  on  the  relative  strengths  of  the  deexcitation 

rates  from  the  3d  excited  states  to  those  of  DR  from  the  fluorine-like  states 

but  also  on  the  relative  ion  populations  that  are  calculated  for  these  states. 

In  CRE,  on  the  other  hand,  we  find  there  is  a  large  abundance  of  ions  in  the 

neon-like  ground  state  and  that  they  contribute  strongly  to  the  excitation  of 

the  JeO  levels.  However,  in  CRE,  we  find  that  the  pumping  of  the  J=2  states  is 

due  predominantly  to  deexcitation  of  the  3d  states,  and  that  it  is  also 

competitive  with  the  collisional  excitation  of  the  J^O  states.  These  points 

are  illustrated  in  Figs.  6  and  7,  which  show  the  calculated  pumping  strengths 

in  CRE  to  the  (1/2, 1/2)  J=0  and  (1/2, 3/2)  J=2  lasing  levels  as  a  function  of 

19  3 

temperature  at  the  fixed  ion  density  of  10  ions/cm  . 

When  we  compare  Figs.  6  and  7,  we  see  that,  even  though  the  contribution 
of  DR  pumping  to  the  J*2  state  is  comparatively  much  larger  than  its  effect  in 
populating  the  J=0  level,  the  strengths  of  DR  pumping  to  both  of  these  levels 
are  much  smaller  than  collisional  excitation  from  the  ground  state.  Fig.  2  of 
Ref.  19  shows  how  dielectronic  recombination  feeds  the  upper  and  lover  lasing 
levels  and  how  the  J=2  levels  can  be  substantially  populated  directly  from  the 
fluorine-like  ground  state  via  DR  when  there  are  enough  F-like  ions.  The  DR 
rate  coefficients  from  An=0  F-like  to  J=0  and  J=2  upper  3p  levels  and  J=1  3s 
lower  lasing  levels  are  such  that  inclusion  of  these  rates  also  tend  to 
increase  the  rate  coefficients  to  the  J=2  levels  compared  to  that  of  the  J=0 
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levels.  However,  when  we  compare  Fig.  2  of  Ref.  19  to  Fig.  3  of  this  work,  it 
is  evident  that  the  dominance  of  collisional  pumping  from  the  neon-like  ground 
state  over  DR  pumping  from  the  fluorine-like  ground  states  that  is  seen  in 
Fig.  7  is  attributable  to  the  preponderance  of  neon-like  ground  states  to 
fluorine-like  ground  states  in  CRE  under  plasma  conditions  relevant  to 
selenium  X-ray  lasing.  When  we  include  detailed  calculations  of  DR  from 
oxygen-like  to  fluorine-like  states,  this  effect  is  increased.  To  further 
understand  the  magnitudes  of  the  curves  shown  in  Figs.  6  and  7,  we  turn  to 
Figs.  8-11. 

Fig.  8  shows  the  collision  strengths  used  in  our  calculations  for  all  the 
multiplet  levels  of  the  3p  state  that  were  obtained  from  Ref.  14.  This  data 
directly  explains  the  differences  in  the  strengths  of  the  J=0  and  J=2 
collisional  excitation  pumping  rates  seen  in  Figs.  6  and  7.  Since  collisional 
excitation  to  both  J«0  and  J«2  levels  is  from  the  same  neon-like  ground  state, 
the  relative  magnitude  of  the  collision  strengths  to  these  two  states 
determines  the  relative  strength  of  the  pumping  rates.  As  can  also  be  seen 
from  Fig.  8,  the  collision  strengths  to  the  other  3p  levels  fall  off  very 
rapidly  with  energy  and  consequently  have  even  smaller  excitation  rates  than 
do  the  J*0  and  J=2  levels. 

The  situation  regarding  deexcitation  pumping  from  the  3d  states  is  more 
complex.  In  this  case,  the  degree  of  metastability  of  the  3d  states  is  an 
important  factor  in  the  3p-3d  coupling.  If  ve  ignore  this  feature  initially 
and  examine  the  coupling  to  the  3p  states  by  summing  over  the  3d  states,  then 
we  obtain  the  results  shown  in  Fig.  9  for  the  3p-3d  collision  strengths  to  the 
JsO  and  J=2  levels  as  obtained  from  Ref.  18.  It  is  seen  from  this  figure  that 
3d  collisional  couplings  to  the  J=2  3p  levels  are  about  an  order  of  magnitude 
higher  than  those  to  the  J*0  3p  levels.  However,  the  pumping  rates  due  to  3d 
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deexcitation  to  the  J>2  levels  in  Figs.  6  and  7  are  about  two  orders  of 

magnitude  larger  than  the  3d  pumping  rates  to  the  J>=0  levels.  In  order  to 

explain  this  behavior  we  have  to  understand  how  each  of  these  3p  (J=0  and  J=2) 

levels  are  coupled  to  each  of  the  3d  levels  and  also  how  each  of  these  3d 

levels  are  populated  in  CRE.  Any  reduction  in  3d  populations  will  give  lower 

deexcitation  pumping  of  the  3p  states. 

A  detailed  examination  of  the  coupling  strengths  between  the  3d  and  3p 

states  shows  the  pattern  depicted  in  Fig.  10;  namely,  each  3d  state,  while 

collisionally  coupled  to  each  3p  state,  is  predominantly  coupled  to  only  one 

as  indicated  in  Fig.  10.  Moreover,  under  the  plasma  conditions  shown  in  Figs. 

6  and  7,  the  3d  states  are  depopulated  predominantly  by  collisional 

deexcitation  rather  than  by  radiative  decay  to  the  3p  states.  The  J-0,  3p 

3  1 

states,  in  particular,  are  strongly  coupled  only  to  the  two  (  and  P^)  3d 
states  that  strongly  radiate  to  the  ground  state;  whereas,  the  J=2  states  are 
strongly  coupled  to  metastable  3d  states.  (The  rate  coefficients  for  the  two 
3d  to  ground  state  radiative  decays  are  also  shown  in  Fig.  10).  Thus  the  3=2 
3p  levels  receive  substantial  population  from  the  3d  metastable  states 
compared  to  the  J*0  3p  levels  because  the  populations  (as  well  as  the  coupling 
strengths)  of  these  metastable  states  are  much  larger  than  the  populations 
(and  coupling  strengths)  of  the  decaying  J=1  states  that  populate  the  J=0  3p 
levels.  To  quantify  these  statements,  the  relative  fractional  populations,  f^= 

O  J  O  J 

f^  “/Zf^  ,  of  the  five  highest  lying  3d  J  levels  are  shown  in  Fig.  11  as  a 

function  of  ion  density  at  an  electron  temperature  of  1  keV.  This  figure  shows 

that  the  fractional  populations  of  the  metastable  J=3  3d  levels  are 

approximately  an  order  of  magnitude  larger  than  the  populations  of  the  J=1  3d 

19  -3 

levels  at  a  plasma  density  of  10  cm  ,  and  that  the  3d  populations  only 

21  -3 

begin  to  equilibriate  at  ion  densities  larger  than  10  cm 
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Ve  see  from  the  above  discussion  hov  the  calculated  populations  of  the 
states  responsible  for  populating  the  laser  levels  play  a  major  role  in 
determining  the  relative  and  absolute  population  inversions.  Thus  it  is 
crucial  to  determine  these  populations  as  accurately  as  possible  for  the  time 
dependent  plasma  conditions  that  are  generated  in  each  experiment.  Accurate 
comparisions  between  theoretically  calculated  and  experimentally  inferred  gain 
coefficients  require  accuracy  in  each  of  these  three  calculations:  rate 
coefficients,  ion  populations,  and  plasma  conditions.  For  this  reason,  X-ray 
lasers  provide  an  important  challenge  to  our  theoretical  ability  to  calculate 
and  comprehend  the  plasma  dynamics  necessary  to  produce  them  efficiently. 
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IV.  GAIN  CALCULATIONS 


It  follows  directly  from  the  discussion  of  the  previous  section  that  in 
order  to  produce  the  maximum  gain  for  the  important  lasing  lines  one  must 
generate  the  right  plasma  conditions  of  temperature,  density,  and  ionization 
where  an  abundance  of  neon-like  or  fluorine-like  ground  state  Se  ions  exist. 
In  CRE,  the  total  fractional  abundances  of  Na-,Ne-,  and  F-like  selenium  ions 
as  a  function  of  electron  temperature  at  an  ion  density  of  10  cm  that  are 
calculated  by  our  model  are  shown  in  Fig.  12.  These  populations  were  obtained 
by  summing  over  all  states  in  each  ionization  stage.  Ve  note  from  this  figure 
that  only  at  electron  temperatures  well  beyond  1  keV,  the  temperature  most 
relevant  to  the  Ne-like  Se  x-ray  laser,  do  the  populations  of  F-like  ions 
become  comparable  to  and  larger  than  those  of  the  Ne-like  ion  in  CRE. 
Nevertheless,  it  is  important  to  calculate  and  include  accurate  atomic  data 
for  the  F-like  ion  from  mechanisms  such  as  excitation,  ionization  and 
recombination  in  order  to  accurately  calculate  the  neon-like  as  well  as  F-like 
populations. 

One  way  to  evaluate  the  quality  of  the  CRE  calculations  resulting  from  our 
model  is  given  in  Fig.  13.  In  this  figure,  the  gain  curves  that  are  computed 
from  the  CRE  population  inversion  generated  between  the  total  number  of  3p  and 
3s  states  in  each  of  the  F-like  and  Ne-like  ionization  stages  are  compared  as 
a  function  of  density  at  an  electron  temperature  of  1  KeV.  To  obtain  the  Ne- 
like  curve,  we  summed  over  the  individually  computed  populations  of  the  3s  and 
3p  states.  In  this  and  all  the  other  gain  curves  to  be  shown  in  this  section, 
we  used  a  Doppler  linewidth  calculated  by  assuming  that  the  ion  temperature 
equalled  the  electron  temperature  T^.  In  Fig.  12  the  3s-3p  gain  coefficient 
for  the  Ne-like  ion  is  much  larger  than  that  for  F-like.  This  result  is 
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and  with 


24 

consistent  with  recent  gain  calculations  of  P.  Hagelstein 

experimental  observations.  Hagelstein  calculated  gain  coefficients  between  2 

and  4  cm~^  for  the  four  3s-3p  J-J  transitions  in  F-like  selenium.  No 

amplification  of  these  lines  has  yet  been  measured,  however,  which  is  in 

keeping  with  Fig.  13,  where  the  calculated  gains  in  F-like  Se  are 

significantly  smaller  than  those  calculated  for  neon-like  Se  at,  for  example, 
19  3 

10  ions/cm  .  Hagelstein' s  calculation  was  performed  for  the  same  plasma 

3 

conditions  as  quoted  in  the  Livermore  experiment  i.e.,  at  the  electron 
20  -3 

density  of  5x10  cm  .  Note,  also,  in  Fig.  12  that  the  gain  in  the  F-like  ion 

peaks  and  falls  off  at  a  slightly  higher  density  than  that  of  the  Ne-like  ion. 

Thus  to  observe  amplification  of  F-like  lines,  one  may  require  plasma 

densities  higher  than  those  yet  achieved.  On  the  other  hand,  while  the  results 

in  Fig.  13  are  suggestive,  they  are  not  definitive.  To  obtain  a  more  accurate 

description  of  the  population  dynamics,  one  must  calculate  the  gain  for  each 

of  the  specific  transitions  among  the  F-like  J  levels. 

Fig.  14  shows  the  gain  coefficients  that  are  calculated  for  the  four  J-0-1 

+24 

and  J«2-l  lasing  transitions  in  Se  as  a  function  of  temperature  at  an  ion 
19  -3 

density  of  1x10  cm  .  The  coefficients  are  labeled  according  to  Fig. 2.  The 

temperature  dependence  of  the  gain  was  calculated  at  this  particular  ion 

density  because  it  is  of  the  same  order  of  magnitude  as  the  density  that  was 

inferred  in  the  Livermore  experiment.  Because  the  gain  coefficients  of  these 

four  J«0-1  and  J*0-2  lines  all  peak  at  1  keV  at  this  density,  we  used  this 

temperature  to  calculate  the  density  dependence  of  these  gain  coefficients, 

which  is  plotted  in  Fig.  15.  Fig.  15  shows  that  the  gain  curves  for  all  the 

19  -3 

four  lines  peak  at  a  larger  ion  density  than  1x10  cm  and  that  the  gains 
for  different  lines  peak  at  different  densities.  Thus  if  all  four  lines  are 
amplified  during  an  experiment,  one  may  infer  an  upper  bound  on  the  density 
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depending  on  the  reliability  of  the  atomic  model. 

While  the  sizes  of  the  gain  coefficients  in  Figs.  14  and  15  are  in  general 

agreement  with  recent  calculations  of  London  et.  al.  at  the  temperature  and 

density  of  their  calculation,  there  is  one  noteworthy  difference.  In  our 

calculation,  only  the  upper  of  the  two  J=0-1  lines  has  a  larger  gain  than  its 

neighboring  J=2-l  line.  Our  calculated  gain  for  the  169  A°  line  is  small  and 

in  general  accord  with  the  early  experimental  observations.  Only  the 

calculated  gain  of  the  182  A°  line  would  appear  to  suggest  problems  in  our 

model  because  of  its  unusally  large  size. 

Because  the  J=0-1  182  A^  line  is  calculated  to  have  a  much  larger  gain 

than  the  J*2-l  lines,  which  may  contradict  experimental  observations,  we  were 

led  to  investigate  the  sensitivity  of  this  gain  coefficient  to  the  strength 

of  the  monopole  collisional  excitation  rate.  Moreover,  there  are  other 

instances  in  the  literature  where  the  accuracy  of  calculated  monopole  rates 

25 

have  been  called  into  question.  It  was  shown  by  Machado  et  al  ,  for  example, 

that,  while  the  calculated  differential  and  integral  cross  sections  for 

electron-impact  excitation  of  all  2p^3s  levels  were  very  close  to  experimental 

26 

measurements,  a  strong  disagreement  with  experimental  data  could  be  found 

for  the  cross  sections  of  the  2p^3p  levels  of  neon.  In  a  different  calculation 

27 

of  the  excitation  of  argon,  Bubelev  et  al  showed  that  the  cross  sections  for 

excitation  of  3p^4p  J*2  levels  had  quite  good  agreement  with  experimental 

28  5 

data  ,  but  that  the  cross  sections  of  the  3p  4p  J=0  levels  were  calculated  to 

have  an  order  of  magnitude  higher  values  than  those  measured  experimentally. 

29 

In  another  example,  Griffin  et.  al  calculated  the  electron  impact  ionization 
cross  sections  of  Na-like  ions  using  distorted  wave  approximation  in  order  to 
investigate  the  effect  of  excitation-autoionization  on  these  cross  sections. 
Their  calculations  indicated  that  the  largest  contributions  came  from  the  2p- 
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3p  excitations  for  all  the  ions  they  considered.  Again,  when  they  compared 
their  results  with  experimental  data,  there  vas  substantial  disagreement.  They 
concluded  that  the  primary  discrepancy  arose  due  to  their  overestimation  of 
the  magnitude  of  the  monopole-dominated  2p-3p  (and  to  a  lesser  extent  the  2p- 
4p)  transition.  All  this  evidence  suggests  that  it  is  plausible  that  the 
calculated  monopole  excitation  rate  coefficient  to  the  upper  J=0  level  that  we 
employed  in  our  model  could  also  be  in  error  by  at  least  a  factor  of  three. 

In  order  to  test  this  hypothesis,  we  reduced  the  collisional  excitation 

rate  to  the  upper  J=0,  3p  level  by  a  factor  of  three  and  used  this  modified 

rate  in  the  CRE  calculation  to  repeat  the  gain  calculations  of  Fig.  15.  As 

expected,  the  gain  for  the  J=0-1  182A®  line  was  reduced,  but,  additionally, 

the  reduction  was  significant.  The  gains  for  the  four  J=0-1  and  J=2-l  lines 

using  the  reduced  monopole  rate  to  the  (1/2, 1/2)^  level  is  shown  in  Fig.  16  as 

a  function  of  density  at  an  electron  temperature  of  1  keV.  Comparing  Fig.  16 

with  Fig.  15,  we  see  that  the  strengths  of  the  gain  coefficients  for  both  sets 

of  J-0  and  J«2  levels  are  now  in  the  same  relation.  The  J=2-l  line  near  206A° 

now  has  the  highest  overall  gain,  and  the  gain  curves  for  the  two  J=0-1  lines 

are  very  close  to  one  another.  The  relative  strength  of  the  gain  curves  in 

this  figure  would  appear  to  be  in  much  closer  accord  with  the  early 

expermental  observations  than  any  of  the  calculations  performed  to  date.  By 

contrast,  if  we  increase  the  collision  excitation  rate  to  the  lower  of  the  two 

J-0  levels  by  only  a  factor  of  two,  the  gain  ratio  for  both  of  the  J=0-1  to 

20  3 

the  J=2-l  lines  at  densities  above  10  ions/cm  can  be  reversed  from  that  in 
Fig.  16.  This  effect  is  shown  in  Fig.  17. 

From  the  results  of  these  gain  calculations,  it  is  apparent  that  one  must 
have  accurate  monopole  collisional  excitation  rates  in  order  to  accurately 
calculate  gain  ratios  involving  the  J=0  levels  of  the  neon-like  3p  states.  It 
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is  also  apparent  that  a  realistic  calculation  of  the  ionization  balance  will 
rely  on  an  accurate  knowledge  of  plasma  conditions.  For  these  reasons,  it  may 
be  easier  and  more  meaningful  to  make  experimental  and  theoretical  comparisons 
between  the  gain  ratios  of  lines  than  it  is  to  compare  the  absolute  values  of 
calculated  and  measured  gains. 


113 


V.  SUMMARY 


The  problem  of  comparing  calculated  with  experimentally  inferred  gain 

coefficients  in  order  to  confirm  our  understanding  of  the  dynamics  of  X-ray 

laser  schemes  is  complex  yet  decipherable  in  general.  The  complexities  are 

due,  in  part,  to  the  fact  that  experimental  plasma  conditions  vary  rapidly 

both  spatially  and  temporally,  and  the  x-ray  laser  beam  may  be  refracted, 

partially  absorbed,  or  receive  non-uniform  amplification.  In  a  recent  pair  of 
30  31 

papers  *  ,  it  vas  suggested,  moreover,  (1)  that  the  observed  ratio  of  gains 

between  the  J«0-1  and  3^2-1  lines  of  neon-like  germanium  would  vary  depending 
on  the  geometry  that  was  used  to  carry  out  the  X-ray  laser  experiment  and  (2) 
that  the  average  temperatures  or  densities  reached  in  X-ray  laser  experiments 
were  not  necessarily  those  predicted  by  hydrodynamics  calculations.  To  this 
list  of  potential  difficulties  in  making  calculational/  experimental 
comparisions,  one  can  add  the  problems  of  (1)  calculating  accurate  rate  data 
for  each  atomic  process,  (2)  building  a  complete  atomic  model  from  which  both 
excited  state  populations  and  ionization  abundances  can  be  self-consistently 
obtained  and  (3)  accurately  calculating  the  ion  populations  for  different 
specific  sets  of  rapidly  changing  plasma  conditions.  In  spite  of  all  these 
difficulties  ,  however,  there  are  several  issues  concerning  the  behavior  of  X- 
ray  laser  gain  coefficients  that  can  be  fruitfully  examined  jointly  by 
experiments  and  by  a  variety  of  supporting  calculations.  One  involves  the 
design  of  experiments  to  promote  one  pumping  mechanism  over  another  in  order 
to  preferentially  amplify  one  set  of  lines  over  another.  Another  involves 
experiments  designed  to  explore  the  density  dependence  of  the  gain 
coefficients  of  different  lines.  Finally,  a  third  involves  experiments 
systematically  carried  out  to  determine  the  atomic  number  scaling  of  the 
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relative  strengths  of  the  gains  of  neighboring  lines.  In  this  paper,  ve  have 
presented  results  of  gain  calculations  in  neon-like  selenium  that  bear  on 
these  latter  two  issues  and  examined  in  detail  the  relative  strengths  of  the 
neon-like  pumping  mechanisms  in  CRE  that  bear  on  the  first  issue. 

A  complete  selenium  atomic  model  consisting  of  30  neon-like  and  eight 

lumped  F-like  excited  states  was  used  to  calculate  self-consistent 

populations,  pumping  rates,  and  gain  coefficients  for  the  3s-3p  lasing 

transitions  in  CRE.  The  model  utilized  state-of-the-art  collision,  DR  and  RX 

rates  in  both  Ne-like  and  F-like  ionization  stages  and  was  able  to  reproduce 

all  of  the  essential  aspects  (with  some  important  differences)  of  previous 

gain  calculations  for  all  of  the  3s-3p  lasing  transitions.  The  collisional 

excitation  rates  to  the  n=3  states  were  the  same  as  those  used  by  the 

Livermore  researchers  and  our  DR  rates  compared  very  well  with  those  of 
20 

M.Chen  (although  we  could  only  compare  total  ground-to-ground  DR  rates  with 
those  of  Chen).  Ve  found  that  the  inclusion  of  state  specific  DR  rates  had 
only  a  minor  direct  impact  in  CRE  on  the  gain  calculation  of  the  3^2-1  lines 
presumably  in  contrast  to  Livermore  findings  .  Even  though  these  rates  were 
calculated  in  a  very  accurate  and  detailed  way,  the  variation  of  the  gain 
curves  as  a  function  of  the  strength  of  the  J=2  3p  DR  rates  was  not 

significant  when  these  rates  were  varied  by  up  to  a  factor  of  three.  The 

reason  for  this  insensitivity  was  that  in  our  CRE  calculation,  collisional 
excitation  and  3d  deexcitation  pumping  dominates  DR  pumping  and  so  tne  latter 
will  have  a  major  influence  only  in  a  non-equilibrium  plasma. 

Another  apparent  difference  between  our  results  and  other  published 
results  was  that  our  model  predicted  the  gain  ratio  J=  (1/2, 1/2)^  -  J= 
(1/2, 1/2)^  /  J=(l/2, 3/2)2  -  J=  (l/2,l/2)j  to  be  larger  than  one  while  J= 
(3/2, 3/2)^  -  J-  (3/2,l/2)j  /  J=(3/2, 3/2)2  -  J=  (3/2, 1/2)^  was  less  than  one. 
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However,  a  decrease  in  the  size  of  the  calculated  (1/2, 1/2)^  monopole 

excitation  cross  section  by  only  a  factor  of  three  was  sufficient  to  reverse 

the  J-  (1/2, 1/2)^  -  J-  (1/2, 1/2)^  /  J-(l/2, 3/2)2  '  (1/2, 1/2)^  ratio  to 

19 

values  less  than  one.  Moreover,  in  this  case,  the  low  density  (<5x10 
.  3 

lons/cm  )  gains  for  both  3=2-1  lines  became  approximately  equal.  The  same  was 

true  for  the  J»0-1  lines.  This  is  an  important  finding  since  it  is  generally 

believed  that  cross  section  calculations  are  no  more  accurate  than  a  factor  of 

two  or  three.  Equally  important,  we  pointed  out  that  this  is  not  the  first 

instance  when  the  size  of  a  calculated  monopole  cross  section  was  suspect.  In 

a  previous  instance,  ’  it  was  believed  that  the  calculated  monopole  cross 

section  was  too  large  by  a  factor  of  ten.  Ve  hope  that  these  findings  will 

provide  sufficient  motivation  for  other  researchers  to  more  accurately 

recalculate  or  measure  the  rates  for  processes  exciting  the  J=0  3p  states.  The 

results  of  our  calculations  would  suggest  that  the  collisional  excitation  rate 

for  populating  the  upper  level  of  the  mysterious  "J=0-1"  182  A°  line  needs  the 

most  attention,  either  theoretically  or  experimentally.  A  thorough 

32 

measurement  of  this  rate  using  the  electron  beam  ion  trap  (EBIT)  at  Lawrence 
Livermore  Laboratory,  for  example,  would  provide  an  important  benchmark  for 
theoretical  atomic  models  of  the  neon-like  system. 

A  second  method  for  benchmarking  gain  calculations  against  experiment  is 
provided  by  the  predicted  density  behavior  of  different  models.  Without  any 
modifications  of  the  monopole  excitation  cross-sections,  the  gain  of  the  J=0-1 

line  at  182  A°  peaked  above  lOOcm”^  at  3x10^^  ions/cm^,  and  it  had  values 

-1  213 

larger  than  10  cm  up  to  densities  of  2x10  ions/cm  ,  well  past  those 

densities  where  the  other  lines  can  be  amplified.  These  results  are  very 

sensitive  to  the  strength  of  the  (1/2,1/2)q  monopole  excitation  rate  since  a 

factor  of  three  reduction  in  its  magnitude  caused  the  peak  gain  of  the  182  A° 
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line  to  be  reduced  by  a  factor  of  5  and  its  gain  to  vanish  above  5x10 

3 

ions/cm  .  Another  important  feature  of  our  calculations  was  that  th  >  gain  of 
the  210  A?  line  vanished  near  or  below  2x10^^  ions/cm^  in  all  of  the 

cases  ve  studied.  This  behavior  suggests  that  if  experiments  could  be  carried 
out  in  which  the  density  of  the  gain  medium  could  be  systematically  increased 
(with  care  to  avoid  refraction  problems),  then  not  only  could  greater 
amplification  be  achieved  but  variations  in  the  relative  strengths  of  the 
amplified  lines  could  also  be  used  to  test  the  accuracy  of  the  rate  data  used 
in  the  calculations. 
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FIG.  1  Energy  level  diagram  for  Ne-like  selenium  as  used  in  our  atomic  model. 


J  =  0-1  AND  J  =  2-1 

X-RAY  LASER  TRANSITIONS  IN  NE-LIKE  SELENIUM 


2p®3s  2p53p 

FIG.  2  Energy  level  diagram  for  the  3s  and  3p  multiplet  levels  of  neon-like 
selenium.  The  important  lasing  lines  are  labeled  and  shown  along  with  their 
corresponding  wavelengths. 
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PIG.  8  Collision  strengths  of  the  ten  3p 
the  ground  state  of  Ne-like  selenium  as  a 
to  the  threshold  energy.  The  circles  are 
(Ref.  14). 
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FIG.  9  Total  collision  strengths  of  the  important  four  excitations  from  the 
J=0  and  J=2  levels  of  the  3p  state  to  the  lumped  3d  state.  The  circles  are 
data  points  from  Hagelstein  and  Jung  (Ref.  18). 
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3d  J=l 


FIG. 10  Energy  diagram  shoving  the  strongest  pumping  channels  of  the  3p 
multiplet  levels  by  deexcitation  from  the  3d  levels.  The  channels  coupling  to 
the  two  strongly  decaying  3d  states  are  indicated  by  lines  with  dark  dots. 
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TOTAL  FRACTIO 


FIG. 12  Total  fractional  populations  of  F-,  Ne-,  and  Na-like  selenium  ions  as 
a  function  of  electron  temperature  at  an  ion  density  of  10^^  cm”^. 
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PIG. 14  Gain  coefficients  (in  cm"  )  for  the  four  J=0-1  and  J=2-l  3s-3p  lasing 
transitions  of  neon-like  selenium  as  a  function  of  electron  temperature.  The 
curves  are  labeled  as  in  Fig.  2. 


T  =  keV 


FIG. 16  Gain  coefficients  (in  cm  )  of  the  J*0-1  and  J=2-l  3s-3p  lasing 
transitions  as  a  function  of  ion  density  at  an  electron  temperature  of  1  keV. 
These  coefficients  were  calculated  using  a  modified  atomic  model  in  which  the 
magnitude  of  the  monopole  excitation  rate  to  the  (1/2, 1/2)^  level  was  reduced 
by  a  factor  of  three. 
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PIG. 17  Gain  coefficients  (in  cm”  )  of  the  J=0-1  and  J=2-l  3s-3p  lasing 
transitions  as  a  function  of  ion  density  at  an  electron  temperature  of  1  keV. 
These  coefficients  were  calculated  using  a  modified  atomic  model  in  which  the 
magnitude  of  the  monopole  excitation  rate  to  the  (3/2, 3/2)^  level  was 
increased  by  a  factor  of  two. 


